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MIDCOURSE -GUIDANCE PROCEDURE WITH SINGLE POSITION FIX 


OBTAINED FROM ONBOARD OPTICAL MEASUREMENTS 

By Harold A. Hamer, Katherine G. Johnson, 
and W. Thomas Blackshear 
Langley Research Center 

SUMMARY 

A manual procedure is developed for midcourse-guidance application. The pro- 
cedure requires only limited onboard calculations and leads to guidance predictions suffi- 
ciently accurate for emergency or backup operations. Although high accuracy cannot be 
obtained by manual methods as compared with automatic methods, the accuracy of the 
present method is shown to be adequate to the point where a terminal guidance system 
could take over with relatively simple equipment. The method is applied to Earth-Moon 
trajectories and is studied in detail for an example lunar mission. 

Basically, the determination of only one position fix is required rather than the 
determination of a large number of position fixes as is usually required in any manual 
or automatic procedure. This condition is made possible by use of precalculated data 
which linearly relate deviations from the nominal trajectory at the time of the position 
fix to those farther along the trajectory. This procedure is especially important for 
lunar missions where time is limited for making measurements and calculations prior to 
the first midcourse maneuver. The effects on aim-point accuracy of measurement 
errors, maneuver errors, star selection, and linear approximation are discussed. 

INTRODUCTION 

The navigation and guidance operations in manned space missions will normally be 
accomplished by elaborate statistical procedures (such as those described in refs. 1 
and 2), which require extensive use of automatic computing equipment to process large 
quantities of onboard and/or Earth-based measurements. A procedure of this kind is 
necessary in order to obtain the accuracy requirements dictated by manned space flight. 
In case of malfunction or failure of this complex system, it would be highly desirable to 
have a backup procedure wherein the astronaut could, from several onboard measure- 
ments, manually guide his spacecraft by the use of relatively simple guidance calcula- 
tions. (Hereafter in this report, navigation and guidance are referred to merely as 
guidance.) 



A method that could be applied to a manual guidance procedure is developed herein. 
The method is based on a fixed-time-of-arrival condition; that is, the spacecraft is 
guided to a certain point on a nominal or reference trajectory. Although the method is 
developed for Earth-Moon trajectories, it could serve for interplanetary trajectories. 

The method pertains to midcourse guidance — that is, guiding the spacecraft in a trans- 
lunar trajectory from a point relatively near the Earth to a point near the Moon. As 
would any manual method involving such large distances, the present method exhibits 
relatively large errors at the aim point (inaccuracies of several hundred kilometers for 
present-day instrument angular measurements) and would normally be intended only as 
an emergency procedure to get the spacecraft back near the right course. Additional 
refinement of the trajectory, which may be necessary farther along the trajectory, could 
be achieved by manual terminal guidance methods such as those described in refer- 
ences 3 to 5. As an example, for a manned circumlunar mission in which the spacecraft 
is to circumnavigate the Moon and return directly to Earth, the present method could 
serve to correct the translunar Earth-Moon trajectory and a terminal method could be 
used to insure a reasonable perilune at the Moon and a safe perigee (or corridor) for 
Earth reentry. 

The method developed herein is based on a linear perturbation procedure, which 
makes use of the so-called transition matrix, and permits much of the computation to be 
handled by preflight calculations. The method is developed from one that makes use of a 
knowledge of vehicle position at two discrete times along the trajectory to calculate the 
velocity and hence to determine the magnitude and direction of the velocity correction 
required for the midcourse maneuver. Although repeated range measurements may be 
required for increased accuracy, the method basically employs a single position fix 
relating to a preselected time. The position fix is obtained from onboard optical angular 
measurements between certain combinations of stars, the Earth, and the Moon. The off- 
nominal vehicle -position deviations at different times can be linearly related; therefore, 
measurements are not required for a position fix at the second time. Hence, this posi- 
tion fix can be empirically determined in terms of the first one and most of the error 
associated with guidance from measurements at two discrete points can be eliminated. 

The accuracy characteristics of the method are thoroughly discussed; results of a 
Monte Carlo analysis of measurement errors and the effects of star selection are 
included. Also discussed are effects on guidance accuracy due to the linear approxima- 
tion made in using transition-matrix theory. The range measurement is the most criti- 
cal from the standpoint of accuracy in any type of guidance system; therefore, the various 
means for manually determining range are examined, along with their relative errors. 

In this report the term "manual” Infers that guidance measurements are made 
onboard the space vehicle and that all necessary guidance calculations are performed 
without resort to an automatic computing machine with built-in programs and procedures. 
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However, it is not meant to preclude the use of a simple machine to perform desk-type 
calculations; in fact, such a machine would be the minimum capability required in any 
manual guidance calculation. 


SYMBOLS 

A angle formed at vehicle by line of sight to Earth center and line of sight to 

Moon center 

a constant of ellipse (appendix B); semimajor axis (appendix C) 

AA,AB,AC,A6 difference between actual value and nominal value of A, B, C, and 6 , 
respectively (for example, A A = Aa - An) 


B 


D,E,F 


angle formed at Earth center by line to vehicle and line to Moon center 
constant of ellipse (appendix B) 

perpendicular distance between asymptote of hyperbolic lunar approach 
trajectory and center of Moon 

angle formed at Moon center by line to vehicle and line to Earth center 

distance between center of ellipse and center of circle (appendix B) 

off-nominal position component in direction of star 1, star 2, and star 3, 
respectively (D in figs. 7 and 8 is component in direction of specific stars 
named) 


6D,6E,6F difference between values of D, E, and F, respectively, at Tpf and Tj 
^for example, 6D = D-pj - ^Tpf^ 

DC,RA angle of declination and of right asaension, respectively 

matrices in guidance equations 

ff focus factor 

[G] orthogonal rotation matrix ^eigenvectors of [Scr]) 
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h 


vector perpendicular to instantaneous Earth-Moon-vehicle plane (fig. 31) 


i,j,k unit vector along X-, Y-, and Z-axis, respectively 

Kj,K 2 ,K 3 constants in equations (9) and (10) 

l,m,n direction cosine of line of sight to star with respect to X-, Y-, and Z-axis, 

respectively 

FpI transformation matrix 




P11P12P13 

^ 21 ^ 22^23 

P31P32P33 

q 


Ar 

s 


s’ 


T 


Tl 

Tf 


u 

u’ 


^ elements of |^pj 

constant in equations (18) and (19) 
range or distance, (x^ + 
incremental range or distance, r^ - rj^ 

position deviation from nominal trajectory, (ax^ + Ay^ + Az^) 

position deviation from trajectory resulting from midcourse maneuver 
performed with no execution errors 

time from injection 

time of position fix 

time of first midcourse guidance maneuver 
time of second or final (aim-point) maneuver 

velocity deviation from nominal trajectory, ^Ax^ + Ay^ + Az^^^^ 

velocity deviation from trajectory resulting from midcourse maneuver per- 
formed with no execution errors 
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V 


velocity 


AV 

AVf 

6V,60',6(^>’ 

X,Y,Z 

x,y,z 

Ax,Ay,Az 


x,y,z 

Ax,Ay,Az 


6x,6y,6z 

X',Y',Z’ 

x',y’,z* 


first-midcourse -maneuver velocity correction 
final- or second-midcourse-maneuver velocity correction 
random errors in equation (A2) 

rectangular right-hand axis system in which X-axis is in the direction of 
Aries, XY-plane is parallel to Earth equatorial plane, and Z-axis is in the 
direction of celestial north pole 

position coordinates in X, Y, Z system 

off-nominal position component in direction of X-, Y-, and Z-axis, respec- 
tively — for example. Ax = Xa - x^ (the notation ^Ax} represents the 



velocity coordinates in X, Y, Z system 


off-nominal velocity component in direction of X-, Y-, and Z-axis, 

respectively — for example. Ax = Xa - Xn notation {^Ax^ represents 



velocity error at injection (appendix A) 
vehicle-centered normal axes of error ellipsoid 
position coordinates in X', Y', Z' system 
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X",Y",Z" rectangular right-hand axis system (appendix B) 


x'’,y",z" position coordinates in X", Y", Z" system 

X'",Y'",Z'" rectangular right-hand axis system (appendix B) 

x'",y’",z"' position coordinates in X"’, Y'", Z’" system 

Xr,Yr,Zr rotating rectangular right-hand axis system in which Xf-axis lies along 

Earth-Moon line, X^Yr-plane is in Earth-Moon plane, and Zj—axis is in 
northerly direction 

xr,yr,Zr position coordinates in Xr, Yy, Zr system 

a,y half-cone angles for star-to-body measurements (appendix B) 

j3 angle (or supplement) formed by lines of sight to two stars (appendix B) 

6 angle formed at vehicle by line of sight to star and its projection in the 

instantaneous Earth-Moon-vehicle plane 

€ eccentricity of orbit (appendix C) 

?,X in-plane and out-of-plane pointing direction error, respectively, of 

midcourse -velocity guidance vector 

6 angle formed at vehicle by line of sight to star and line of sight to body 

center (referred to as star-to-body angle); true anomaly (appendix C) 

0j,0jj angular measurements used for determining d (fig. 15) 

6oo angle defined in figure 38 

orientation angles (appendix A) 

0 angle formed at vehicle by line to Earth and projection of line to star in the 

instantaneous Earth-Moon-vehicle plane 

ji product of universal gravitational constant and mass of Moon (appendix C) 
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[^r] 

[S«r'] 


■ cr 




covariance matrix of position errors in X, Y, Z system 
covariance matrix of position errprs in X', Y’, Z' system 
standard-deviation value or 1 -sigma error 
state -transition matrices 


3X3 submatrices in state -transition matrix for transitioning from 
Ti to Tf 

l»r^2l’r^3l>r‘^4l 3x3 submatrtces in state -transition matrix for transitioning from 

L JL JL JL J 


Subscripts: 


A angle formed at vehicle by line of sight to Earth center and line of sight to 

Moon center 

a actual or measured (assuming no error) value 

B angle formed at Earth center by line to vehicle and line to Moon center 

calc calculated from equation (14) 

D,E,F off-nominal position component in direction of star 1, star 2, and star 3, 
respectively 

em Earth center to Moon center 

ev Earth center to vehicle 

n nominal value 

p perilune 

r range measurement 

s position deviation; sphere of influence (appendix C) 
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Tpf 

Tl 

Tf 

u 

V 

V 

ve 

vm 

x,y,z 

x',y’,z' 

x,y,z 

0 % 0 > 

1,2,3 

Notation: 

[] 

T 

[] 

[f 

0 
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at time of position fix 
at time of first midcourse maneuver 
at time of second or final maneuver 
velocity deviation 
injection velocity 

first-midcourse-maneuver velocity 
vehicle to Earth center 
vehicle to Moon center 

position component in direction of X-, Y-, and Z-axis, respectively 

position component in direction of X'-, Y'-, and Z'-axis, respectively 

velocity component in direction of X-, Y-, and Z-axis, respectively 

measurement of angle formed at vehicle by line of sight to star 1, star 2, 
and star 3, respectively, and line of sight to body center 

orientation angles 

star 1, star 2, and star 3, respectively 

square matrix 
transpose of square matrix 

inverse of square matrix 

column matrix 


I nil I 


nil 


III lll■n■llllll^l ~ ni rmum 
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L J row matrix 

j j absolute value 

A bar over a symbol indicates a vector. 


BASIC METHOD 
Synopsis 

In general, any guidance correction is determined from information on the vehicle 
velocity. In the simplest method, velocity can be determined from vehicle -position mea- 
surements at two discrete times along the trajectory. Such a method, however, would be 
unacceptable because small errors in the position measurements would produce exces- 
sive error in the velocity determination and, hence, in the guidance calculation. 

In an attempt to circumvent this deficiency, the characteristics of a large number 
of perturbed trajectories were analyzed with regard to the position deviations from a 
nominal trajectory. It was found, for a wide range of perturbed trajectories, that the 
deviations in the directions of certain stars were linearly related over a given time 
interval. Use of these stars for the guidance equations permits the position at the sec- 
ond time to be predicted on the basis of the initial position measurement and eliminates 
much of the effect of measurement error. 

The procedure for calculating the midcourse -maneuver velocity correction is 
developed from standard guidance equations. These equations incorporate three star-to- 
body measurements along with a range measurement. The range measurement also nor- 
mally includes one or two star-to-body measurements, with limits as to the location of 
the stars. Various criteria are developed herein for preselecting (before the mission) 
suitable stars for the range measurement and for the guidance equations. 

General Considerations 

Transition matrices .- The present method is based on a linear perturbation pro- 
cedure using the so-called state -transition matrix. This matrix can be used to relate 
the deviations in the state of the trajectory at one time to those at some other time. For 
example. 
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where the state -transition matrix maps the position and velocity deviations (from the 
nominal values) at time Tpf to those at time Ti. The state -transition matrix 


8Xt^ 

8Xtj 

8Xtj 

8Xtj 

8xt^ 

8Xtj 

®^Tpf 


9Zrp 

Tpf 




9yTj 


ayTj 



ay^^ 

axrp 

ipf 




"^Tpt 

^ZTpf 

8ZTj 

SzTj 

8ztj 



8zTj 

8xt . 
^pf 


8zt . 
^pf 


"j'Tpt 


8Xtj 


8kTj 

3xtj 

ai^j 

8Xtj 

8xt . 
Apf 


8zt f 
^pf 

8xt , 
^pf 


8Zm 

Tpf 

9YTj 

Sj-Tj 

®yTj 

9YTi 


9yTj 

9Xrp 

^Pf 


8Zrr, 

•Tpf 

8Xrp 

Tpf 


8zrr 

Tpf 

8zt^ 

a^Tj 

8ztj 

8ztj 

8ztj 

8ztj 

9Xrn 

Tpf 

®yT - 
pf 

OZrp 

^Pf 

9Xrp 

^pf 

®^T f 

pf 

8Zrp 

Tpf 


is strictly applicable only when the equations of vehicle motion are linear. Because of 
the actual nonlinearity of these equations, this procedure is limited to perturbed trajec- 
tories that are reasonably close to the nominal. The transition matrix can be computed 
by several methods (see, for example, refs. 6 and 7) and there exist a number of com- 
puter programs (for example, orbit-determination programs) which have the capability 
of generating this matrix for any time interval along a trajectory. 

Measurement of position .- As in any onboard guidance system, the guidance calcu- 
lations are based on measurements of the position of the spacecraft at given times along 
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the trajectory. These measurements are used to determine the spacecraft velocity 
which, in turn, is used to derive the velocity correction required for a midcourse maneu- 
ver. The combinations of onboard optical measurements chosen for the present method 
are Illustrated in figure 1 . 


Z 



Figure L- Required measurements for position fix. 


Of the many combinations of possible onboard measurements that will give a posi- 
tion fix (ref. 8), the method using three star-to-Earth angular measurements together 
with a range measurement (see fig. 1) is the most applicable to manual position fixing. 
(See ref. 9 for a study of this method as applied to a manual procedure.) Because the 
measurements are made onboard, the range is necessarily determined from angular 
measurements, as is discussed in a subsequent section. 

From the aforementioned measurements, the set of equations (for a vehicle- 
centered coordinate system) which leads to the solution of the position vector of the 
Earth for a given time is 



limini 

I2m2n2 

Ismang 


-1 


fr cos 0i| 
r cos 02^ 
|r cos 03! 


or, for the position deviation of the vehicle off the nominal trajectory. 
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limini 

I2m2n2 

Ismsns 



COS 

cos 

cos 


^l,a 
^2, a 
^3, a 


Tji COS 

cos 

cos 


^l,n1 

®2,n> 

^3,nJ 


By letting 


pa cos 01 
^ra cos 02, a 
pa cos 03 , a 


Tji COS 

cos 

cos 


02, n> 

03, nj 



the set of equations may be written as 



limini 

I2m2n2 

I3m3n3 



( 1 ) 


(2) 


where 1, m, and n are the known direction cosines of the respective stars and D, E, 
and F are the position-deviation components in the direction of the stars as illustrated 
in figure 2. 



Figure 2.- Example of variables involved in determining off-nominal position components. 


Development of Guidance Procedure 

The present method is composed of two guidance maneuvers: (1) the primary 
maneuver, which is executed shortly after Earth injection, and (2) the one made at the 
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aim point, necessitated by the first maneuver. The aim point is defined as that point on 
the nominal trajectory to which the first midcourse maneuver attempts to correct. In 
this report, aim points at two different times are investigated: one at perilune and one 
at a point near the lunar sphere of influence. The purpose of investigating the different 
aim points is explained in a subsequent section. 

First-midcourse -man euver guidance equations .- The equations for the first mid- 
course maneuver are developed according to figure 3. 



Figure 3.- Schematic sketch of guidance procedure. 


The deviations from the nominal trajectory at the final (aim-point) time Tf after 
a midcourse guidance correction (see ref. 6, for example) are 



(3) 


where the transition matrix is the mapping matrix from Tj to the final time Tf and 
the prime denotes the velocity deviation immediately after the instantaneous (practically) 
guidance maneuver. For a fixed -time -of -arrival guidance law, the objective is to arrive 
at the aim point on the nominal trajectory, without regard for the final velocity deviation 


Ax. 


Hence 






where the transition matrix is the mapping matrix from Ti to Tpf. Then, 
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From equation (2) 




-1 


limini 

l2i»2“2 

Ismans 


-1 

. — 




\ 

E 

1 

F 


Tpf 


T-ir 






limini 

12 m 2H2 

Ismsng 





( 6 ) 


Therefore, the required first-midcourse-maneuver velocity-correction vector AV, 
obtained by subtracting equation (5) from equation (6), is 



where [^ej and |^gj are 3x3 matrices which can be precalculated from a knowledge of the 

transition matrices and the direction cosines of the pertinent stars. Since the measure- 
ments D, E, and F pertain to a vehicle reference system, the correction vector AV 
is added to the vehicle velocity vector. 

Second-midcourse-maneuve r guidance equations .- In deriving the equations for the 
first midcourse maneuver, the corrected trajectory is designed to pass through the aim 
point; no attempt is made to control the velocity vector at this point. Hence, the second 
midcourse maneuver merely corrects the velocity error introduced at the aim point by 
the derivation of the first midcourse velocity. If the nominal perilune is the aim point, 
the second-midcourse-maneuver velocity correction could be readily combined with the 
lunar deboost maneuver, if this type of mission were contemplated. 

From equation (3), then, the required second-midcourse -maneuver velocity- 
correction vector is 



or, from equation (4), 


AVf = 


#3 fAxrpjj - 




$4 


#2 


-1 


#1 


(axtiJ 


and, from equation (2), 


AVf = 


#3 


#4j $2j ^1 


— — 

-1 

r~\ 

Ijm^nf 


D 

I2m2n2 

< 


^3“3“3 


[fJ, 


By use of the inversion property of the transition matrix (ref. 6) which is 

- 'N" -H" 

l-H" 

it can be shown that the expression for AVf reduces to 




limini 

-1 


AVf = - 

r 

1^2 J 

I2m2n2 

< 

E 1 


1 

l3“3“3 


F r 
^ J • 


Tl 



( 8 ) 


where j^hj is a 3 x 3 matrix that can be precalculated. As in the case for the first mid- 
course maneuver, the correction vector AVf is added to the vehicle velocity vector. 


Method of Determining Position at Tf 

General con sid erations. - As shown by equation (6), the first-midcourse -maneuver 
velocity-correction vector AV is determined from a knowledge of spacecraft position 
at two times. For this type of method, which uses no statistical procedure to reduce a 
large number of position measurements in determining the guidance velocity, extreme 
measurement accuracy is necessary. For example, if a measurement is made at each 
time, the position-determination errors are equally as likely to be opposite in sign 
(direction). Even if the position-determination errors are small at the two times, a 
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large error could be introduced (as illustrated in fig. 4) in the calculation of the velocity 
vector of the actual trajectory and thus a large error could result at the aim point. 





Greater accuracy in determining the velocity vector of the actual trajectory could be 
achieved if the position-determination errors at both times were in the same direction 
and also of the same magnitude, as illustrated in figure 5. 


Actual trajectory 

Measured trajectory 

Figure 5.- Sketch showing position-determination errors that are in the same direction at Tpf and Tj. 

Empirical procedure . - In order to essentially satisfy the conditions illustrated in 
figure 5, a procedure is herein developed whereby the position at time is not meas- 

ured but is predicted from the position at Tpj. This prediction is possible since, for a 
number of perturbed trajectories, a change in deviation in the direction of a particular 
star over the interval Tpf to Tj is a linear function of the deviation at time Tpf. 
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Shown in table I is a comparison of results obtained by the present method and by the 
method involving separate measurements at Tpf and Tj. Values shown were deter- 
mined from a Monte Carlo error analysis (presented in a subsequent section) and corre- 
spond to 1-sigma measurement errors of 10 kilometers in the range and 10 arc-seconds 
in the star-to-body angles. 


TABLE I.- COMPARISON OF ACCURACIES OF PRESENT METHOD AND 
METHOD INVOLVING POSITION MEASUREMENTS AT TWO TIMES 



Position error Og 

Velocity error Ou 

Method 

at Tf = 56 hr. 

at Tf = 56 hr. 


km 

m/sec 

Present 

162 

1.2 

Two-position fix 

2535 

17.5 


The characteristics of the nominal trajectory used throughout this report are given 
in table H. This typical 3 -day translunar trajectory is presented in figure 6 in the 
rotating-axis system which shows the relative positions of the Earth, Moon, and vehicle 
at any given time. 


TABLE n.- CHARACTERISTICS OF NOMINAL TRAJECTORY 


Injection position: 

X, km 

y, km 

z, km 

Injection velocity: 

X, km/sec 

y, km/sec 

k, km/sec 

Coordinate system 

Injection date 

Injection Greenwich mean time, hr:min:sec . . 

Julian date 

Perilune distance, km 

Perilune velocity (selenocentric), km/sec . . 

Approximate time to perilune 

Approximate time to lunar sphere of influence 


-2104.8415 

5395.4303 

2874.9062 


-10.519851 

-2.6309587 

-2.0283715 

Cartesian equatorial of date 

March 18, 1968 

22:27:15 

2439934.43559025 

3403.62 

1.98051 

70 hr, 37 min 

56 hr 
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Figure 6.- Nominal trajectory in a rotating-axis system. (T denotes hours from injection.) 


For the present procedure, a number of perturbed trajectories (see table ni) are 
generated from which values of D are calculated at times Tpf and T^. (In the 
guidance equations, D, E, and F represent different stars.) The values of D can 
be calculated from equation (1) or from the equation 

D = l(xa - xn) + m(ya - Yn) + n(za “ zn) 


where 1, m, and n are the direction cosines of the particular star and x, y, and z 
are vehicle-centered coordinates. Plots of 6D as a function of D for given stars 
are obtained (for example, see figs. 7 and 8). For a given star the variation of 6D with 
D provides the relationship 

6D ^ 

°Tpt+Di5'^Tpi = DTj 

or 

KiDTpf = Dt^ 
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TABLE m.- PERTURBED TRAJECTORIES USED IN STUDY 


Perturbed- 

trajectory 

no. 

Injection-position 
oerturbations, km 

Injection-velocity 
perturbations, m/sec 

Ar 

Off-nominal conditions at aim point 
for perfect measurements and 
perfect guidance maneuver 
at first midcourse when — 







Tpf = 9.5 hr)j 
km 

Tf = 

56 hr 

Tf = 70 hr, 37 min 
(nominal perilune) 


Ax 

Ay 

Az 

Ax 

Ay 

Az 


s, km 

u, m/sec 

s, km 

u, m/sec 

1 

-5 






135.87 

3.263 

0.015 



2 

-10 






271.85 

5.925 

.051 

11.387 

1.285 

3 


-5 





-342.54 

2.668 

.021 



4 


-10 





-686.23 

15,661 

.100 

45.028 

10.676 

5 



5 




182.60 

3.218 

.026 
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10 




365.09 

6.020 

.052 

10.178 

1.828 

7 

2.9 

2.9 

2.9 




225.31 

3.394 

> .027 

7,501 

1.544 

8 




9.15 



-840.68 

13.550 

.094 

40.700 

9.589 

9 




-9.15 



836.58 

16.512 

.071 



10 




4.57 



-419.37 

5.045 

.041 

13.804 

3.005 

11 





9.15 


-207.24 

2.413 

0.015 

6.821 

1.696 

12 






1 

o 

175.67 

3.364 

.024 
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10 

-175.00 

1.156 

.008 

3.418 

.843 

14 






-5 

87.75 

1.587 

.011 
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5 

-87.58 

6.124 

.034 
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-7 

-7 

281.58 

2.574 

.025 

7.991 

1.715 

17 





7 

7 

-281.29 

1.557 

.012 



18 





-7 

7 

36,38 

.372 

.003 



19 





7 

-7 

-35.55 

2.025 

.014 
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5.27 

5.27 

5.27 

-696.20 

13.539 

.084 

39.100 

9.621 

21 




4.57 

5.27 

5.27 

-631.79 

11.729 

0.073 



22 





9.15 

5.27 

-299.68 

2.476 

.018 



23 


-6 


6 



-963.57 

17.568 

,116 

52.050 

12.483 

24 

-5 



-10 



1049.16 

22.760 

.188 

23.800 

5.691 

25 

5 



10 



-1055.57 

21.661 

.147 



26 

-10 



-10 



1184.33 

26.002 

.211 

29.669 

7.922 

27 

10 



10 



-1192.06 

29.210 

.197 



28 


-5 




-10 

-166.44 

3.100 

.027 



29 


-10 




-10 

-509,71 

14.412 

.096 

39.704 

9.070 

30 


10 




10 

507.40 

9.625 

.082 



31 



-5 

10 



-1102.85 

22.291 

0.143 



32 



5 

-10 



1095.55 

23.855 

.198 

24,974 

5.840 

33 



-10 

10 



-1286.84 

27.882 

.177 

82.143 

20.390 

34 



10 

-10 



1276.88 

31.989 

.262 



35 



-5 


10 


-409.47 

6.661 

,041 



36 



-10 


10 


-592.59 

12.935 

.078 

36.697 

9.009 

37 



10 


-10 


591.47 

10.507 

.088 



38 


-5 


10 



-1263.70 

27.648 

.175 



39 


-10 


10 



-1609.62 

42,745 

.266 

123.298 

31.174 

40 


10 


-10 



1591.37 

48.743 

,367 




20 




Values of constants Ki, K2, and K3 for three suitable stars are then used to deter- 
mine off-nominal position components D, E, and F at time Tj for any given per- 
turbed trajectory. This method insures that the position-determination errors are 
about the same at times Ti and Tpf. Substituting the constants into guidance equa- 
tions (7) and (8) yields 



Ckid' 


AVf = [h]< 


K2E > 



( 10 ) 


Star Selection 

This section deals with the problem of selecting stars for obtaining linear plots of 
6D as a function of D. Range -measurement accuracy also depends on the selection of 
stars, as discussed in a subsequent section. 

Geometry considerations .- The overall guidance accuracy of the present procedure 
depends on the three preselected stars used for the position measurement at time Tpf. 
Suitable stars give good linear plots of 6D as a function of D and yield a minimum 
amount of scatter. Scatter is equivalent to position error, and the less scatter, the less 
chance of introducing error in calculating the guidance velocity. 

Figures 7(a) to 7(f) are for suitable stars; stars considered in figures 7(g) and 
7(h) are unsuitable. Figures 7(b), 7(e), and 7(f) are for fictitious stars and illustrate the 
linearity characteristics of stars in the particular directions indicated. 

Experience has shown that the greatest degree of linearity is obtained when the 
star yields maximum values of both D and 6D. This fact leads to consideration of 
probable errors in the injection conditions and the resulting error ellipsoid at Tpf. 

The major axis of this error ellipsoid is the most probable locus of deviations from the 
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(a) Rasalhague. 

Figure 7.- Change in off-nominal position component in the direction of the various indicated stars plotted for different perturbed trajectories 
Data apply to nominal trajectory described in table II. (Signs of 0 and 60 are the same except for data plotted below the zero-line) 













iD'l at Tpf, km 


(f) Fictitious star in direction of major axis of error eiiipsoid, 
Figure 7.- Continued. 
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(h) Regulus. 


7.- Concluded. 







nominal trajectory at time Tpf. As illustrated in figure 9, maximum values of D will 
most likely be measured if the chosen star falls on or near this major axis. Maximum 
values for 5D are obtained by measuring stars in the direction of the velocity vector 
(or along the range vector rye which lies close to the velocity vector in the midcourse 
portion of the trajectory). As shown by figure 9, the nominal and perturbed trajectories 
are approximately parallel; hence, measuring stars perpendicular to the trajectory 
(velocity vector) gives no indication of a change in D between Tpf and Ti. 



Figure 9.- Sketch showing geometry for star seiection. (Ail lines shown lie near the Earth-Moon plane.) 

It appears, then, that for the present guidance method, stars should be selected in 
the general region of the cone whose apex is formed by the error-ellipsoid major axis 
and the range vector fye (see fig. 10). 



Figure 10.- Region of stars suitable for guidance measurements. 
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The most suitable stars would lie in the middle of this region since use of these stars is 
the best compromise for obtaining large values of both D and 6D. 

By referring to the data in figures 7 and 8 along with the angles listed in these fig- 
ures for the respective stars, the differences between the plots can be readily explained. 
The stars cited in figures 7(a) to 7(f) all lie within the region described in figure 10 and 
can be considered suitable. For figure 7(e) with the star in the direction of the vehicle- 
Earth line, the actual measurement, of course, would have to be to a star in the opposite 
direction. A star in the direction of the Earth-Moon line (fig. 7(b)) is shown to be about 
the most suitable since this star is almost directly between rye and the major axis of 
the error ellipsoid. Rigil Kentaurus (fig. 7(c)) lies between these two lines also, but is 
inclined about 39*^ to the Earth-Moon-vehicle plane. The stars Rasalhague (fig. 7(a)) and 
Capella (fig. 7(d)) are in the general direction of the vehicle-Earth line. Regulus 
(fig. 7(h)) is nearly perpendicular to both the major axis of the error ellipsoid and the 
range vector fyg (see fig. 11); hence, the values of both D and 6D are small and 
Regulus could not be used for the present nominal trajectory (table n). 



Procyon (fig. 7(g)) is nearly perpendicular to the major axis of the error ellipsoid; hence, 
the values of D are small and Procyon could not be used for the present nominal tra- 
jectory. The spatial direction of the major axis, however, is dependent upon the partic- 
ular nominal trajectory (injection date, injection velocity, etc.) and, as is shown in fig- 
ure 8, Procyon would be a highly desirable star for the nominal 68-hour translunar tra- 
jectory in reference 10. 
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Calculation of ma jor axis of error ellipsoid.- As indicated in the previous discus- 
sion, a knowledge of the direction of the major axis of the position-error ellipsoid is 
important since it serves to define a limit on the star-to-Earth angle measurements. 
The position-error ellipsoid at Tpf can be determined from 


[ser*] = [G][ser][G]'^ (11) 

where is a diagonal matrix whose elements are vari- 

ances in the values of x’, y', and z’, respectively (for example, see ref. 11). These 
variances are a measure of the three normal axes of the error ellipsoid and are easily 
calculated, as they are the eigenvalues of [ser]> the covariance matrix of position errors 

|at Tpf^ in the X, Y, Z system. The matrix joj is the rotational transformation 
matrix such that 


X’ 

y’ 





Tpf 


)Tpf 


( 12 ) 


This matrix is also easily calculated; the normalized eigenvectors of the matrix [SerJ 

form the rows of [g] which are the direction cosines (with respect to the X-, Y-, and 
Z-axes) of the three normal axes of the error ellipsoid. Hence, the major axis of the 
ellipsoid is determined by the largest eigenvalue, and the direction cosines of the major 
axis are determined by the corresponding row (eigenvectors) of the matrix |^gJ. 

Actually the error ellipsoid is six dimensional, if the velocity errors are considered 
also. However, the position-error characteristics of this ellipsoid, as determined from 
equations (11) and (12) used with six-dimensional matrices, are almost exactly the same 
as those determined from the 3x3 position-error covariance matrix. 


The position-error ellipsoid for the nominal trajectory in the present paper 
(table II) was determined at Tpf = 9.5 hours and the results for various injection errors 
are shown in table IV. The program in reference 7, which uses equations similar to 
equation (11), was employed to obtain the covariance matrix of position errors at 
Tpf = 9.5 hours by mapping forward the covariance matrix of injection errors. Because 
the injection errors are dependent upon the type of vehicle used, arbitrary injection-error 
values had to be assumed. Given in the first column of table IV are the various types of 
assumed errors, some of which are defined in appendix A. Except for those types defined 
in appendix A, the 1-sigma values of the magnitude of the position and velocity errors 
were assumed to be 10 kilometers and 10 m/sec, respectively. The elements of the 
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TABLE IV.- ERROR-ELLIPSOID MAGNITUDE AND DIRECTION AT Tpf FOR VARIOUS INJECTION ERRORS 


Types of errors 


Elements of diagonal covariance matrix 
of injection errors 

1 -sigma values for ellipsoid of 
trajectory-state deviations 

Direction cosines of major axis 
of error ellipsoid 

Angle* 

between major 
axis and 7yg, 
deg 

assumed at injection 

km2 

ay2, 

km2 

0 ^, 

km2 

(km/sec) 2 

(km/sec) 2 

(km/seo)2 

Major axis, 
km 

Mean axis, 
km 

Minor axis, 
km 

1 

m 

n 

Spherical position; 
spherical velocity 

33.33 

33.33 

33.33 

0.00003333 

0.00003333 

0.00003333 

930.19 

124.09 

96.69 

0.73604402 

0.57414687 

0.35860079 

39.037 

Spherical position; 
no velocity 

33.33 

33.33 

33.33 

0 

0 

0 

615.05 

92.67 

78.20 

0.74941316 

0.56076837 

0.35202090 

40.184 . 

No position; 

spherical velocity 

0 

0 

0 

0.00003333 

0.00003333 

0.00003333 

698.58 

95.04 

27.72 

-0.72545294 

-0.58438199 

-0.36361478 

38.145 

No position; 
velocity in direction 
of velocity vector 

0 

0 

0 

0.00009094 

0.00000566 

0.00000339 

1132.54 

37.99 

9.28 

0.74264802 

0.56867225 

0.35367463 

39.592 

Spherical position; 
velocity in direction 
of velocity vector 

33.33 

33.33 

33.33 

0.00009094 

0,00000566 

0.00000339 

1288.76 

93,20 

87.05 

0.74417814 

0.56689993 

0.35330345 

39.725 

No position; 

velocity according 
to appendix A 

0 

0 

0 

0.00009948 

0.00009355 

0.00009653 

1205,79 

159.96 

46.99 

-0.72634358 

-0.58352027 

-0.36322046 

38.220 

Spherical position; 
velocity according 
to appendix A 

33.33 

33.33 

33.33 

0,00009948 

0.00009355 

0.00009653 

1353.47 

179.02 

103.91 

0.73118129 

0.57887490 

0.36094007 

38.625 

Position and velocity 
according to 
appendix A 

87.738 

94.827 

106.236 

0.00009948 

0.00009355 

0.00009653 

1593.09 

206.28 

169.83 

-0.73465308 

-0.57402232 

-0.36163959 

38.930 

Position according to 
■ appendix A; 
no velocity 

87.738 

94.827 
1 

106.236 
J 

0 

J i 

1 

0 

0 

1041.41 

163.26 

128.03 

j 

0.74557412 

0.56113183 

0.35951393 

39.882 


^Either angle or supplement of angle given, whichever is smaller. 


resulting covariance matrix (assumed to be diagonal) of the injection errors are also 
shown in this table. 

It is important to note that for the nominal trajectory, the value of Tpf, and the 
injection-error covariance matrices used in this report, table IV shows the direction of 
the major axis of the error ellipsoid to be essentially unaffected by the choice of injec- 
tion errors. The angle between the major axis and the range vector ryg (see fig. 10) 
is approximately 39°. 


RANGE DETERMINATION 

The accuracy of the range measurement is the predominant factor affecting the 
accuracy of the position fix (see ref. 9) and hence the accuracy of the guidance method. 
Because of the importance of this measurement, the various aspects of the onboard 
optical (manual) determination of range are discussed in this section. 

Methods of Optically Measuring Range 

The two basic methods of determining range to the Earth center from onboard 
optical measurements are: (1) the method involving measurement of the angular diam- 
eter of the Earth, and (2) the method involving measurement of two or more angles 
between certain combinations of stars, the Earth, and the Moon. Large errors are inher- 
ent in angtilar diameter measurements, so the former method would ordinarily not be 
considered. 

The latter method has many alternatives, although basically it involves using the 
sine law to solve the triangle shown in figure 12 for the magnitude of fye- 


To star /—Vehicle 



Figure 12.- Variables involved in determining range from vehicle to Earth center. 
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For any given time , rem is known; the angle A can be measured directly. In order 
to solve for |fve| then, only the angle C remains to be determined. Once the required 
angular measurements have been made, the sine-law equation 


_ sin C 
ve - sin A 


^'em 


(13) 


could be solved directly. This procedure involves use of trigonometric tables and use of 
a large number of significant figures, although calculations show that five-place trigono- 
metric functions provide range -determination accuracy within ±4 kilometers. If a desk 
calculating machine is available, these requirements present no problem as to calcula- 
tion time or calculation accuracy. 

Another way of determining range, which lends itself more readily to manual cal- 
culations, pertains to the perturbation procedure (ref. 9) about a nominal trajectory. For 
this procedure the difference in range from the nominal value is 

Ar„e = - AA + rem(cos B cot A - sin B)AB (14) 

sin A 


or in terms of the angle C 


Aryg = -ryg(cot A)AA + ryg(cot C)AC (15) 

All values in equations (14) and (15) are based on the nominal trajectory and AA, AB, 
and AC are the measured deviations from the respective nominal angles. 

Measurement of Angles B and C 

The solution for range ryg by the method involving the measurement of two or 
more angles requires that either B or C be determined. There are essentially 
three ways by which angle B (or C) can be determined from onboard measurements: 

(1) by measuring the declination and right ascension of the Earth (or Moon), (2) by mea- 
suring one star-to-body angle, and (3) by measuring two star-to-body angles. The first 
method has some disadvantages in that the measurements must be referenced to a stable 
platform. Also, a fair amount of calculation is involved in computing the angle. First, 
the direction of the Earth (for determining B) is determined from 
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( 16 ) 


Xy’g 

= cos DC cos RA 

^ve 


= cos DC sin RA > 
rye [ 


Zye 


tye 


= sin DC 




where DC and RA are the angles of declination and right ascension, respectively. 
Then, from the dot product of the vehicle-Earth vector and the Earth-Moon vector, the 
following equation is obtained; 


cos B = 



cos DC cos RA + Ygjjj cos DC sin RA + 2 .^^ 

^em 


sin DC) 


(17) 


The calculations in equations (16) and (17) could be performed before the mission and 
presented in a form such as that shown in figure 13. (Fig. 13 is discussed subsequently 
in connection with the two-star method.) 

If the one-star method were used, |AB| would be approximately equal to |A0( for a 
star near the Earth-Moon-vehicle plane (see fig. 12) and equation (14) could be written as 


Arye - 


2^ 

7”^ AA + qrem(cos B cot A - sin B)A0 


(18) 


where 0 is the angle between the line of sight to the star and the line of sight to the 
Earth center. Likewise, equation (15) could be written as 


Arye = -rye(cot A)AA + qrye(cot C)A0 (19) 

where 6 is the angle between the line of sight to the star and the line of sight to the 
Moon center. The term q in these equations is either ±1, depending on the location of 
the star (see ref. 9). The one-star method, however, is highly restrictive because of 
accuracy considerations. For reasonable accuracy the star must be in or very close to 
the nominal Earth-Moon-vehicle plane, as is shown subsequently. 

In using the two-star method, the angle B (or C) can be measured directly, as 
shown in appendix B. The calculations are quite extensive, but can be done before the 
flight and presented in the form shown in figure 13 for use -in the perturbation procedure 
given by equations (14) and (15). (The data could be plotted in terms of the total angle C 
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Figure 13.- Difference in angle B from nominal value, as obtained from two onboard star-to-Earth measurements at 
^pf ” hours. (Data apply to nominal trajectory in table II and can be used to determine range by the perturbation 

method.) 
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for direct use in eq. (13)). Although the stars are not required to be in a particular 
plane, they are limited (because of the method used to solve the cone equations) by the 
requirement 


i^+y) < 90° 

Thus, the angle (or supplement) between the lines to the two stars must be less than 
(90° - y). This limitation is not very severe since the same star-to-body angles used 
for the range measurement would ordinarily be used for two of the three position- 
component measurements (D, E, and F). As shown in table V and figure 10, the suit- 
able stars for position measurement are necessarily limited to a small region. 

TABLE V.- ANGLES^ BETWEEN VARIOUS COMBINATIONS OF SUITABLE STARS 


star 

Star ^ 

Fictitious 
(Moon-Earth line) 

Capella 

Rigil 

Kentaurus 

Rasalhague 

Fictitious (vehicle -Earth line) 

20.957° 

16.961° 

41.653° 

13.731° 

Fictitious (Moon-Earth line) 


27.850° 

39.081° 

19.859° 

Capella 

27.850° 


27.073° 

30.301° 

Rigil Kentaurus 

39.081° 

27.073° 


52.520° 


^Either angle or supplement of angle, whichever is smaller. 


Accuracy of Range-Measurement Methods 

A summary of the relative accuracy characteristics is given in figure 14 for the 
various methods of onboard optical range measurement. The standard deviation of range 
error o’ryg calculated from equations (14) and (15). For example. 




rem(cos B cot A - sin B)og 


( 20 ) 


where the angular measurements A and B are considered to have random uncorre- 
lated errors, each with a standard deviation of 10 arc-seconds (for the one-star mea- 
surement method). For the method using the two-star measurements, the standard 
deviation for each of the star measurements was assumed to be 10 arc-seconds. Using 
these values in connection with figure 13 gives an overall 1-sigma error for the angle B 
of 
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(j-Q = ^(14)2 + (14)^ = 19.8 arc-seconds 

The two stars used in figure 13 were Rasalhague and Rigil Kentaurus; however, the 
accuracy shown is representative of any two stars. The reading error is not included 
since it would be relatively small and have little overall effect on the value of Og. It 
should be stated that, in reference to equation (20), error in the angle A has a com- 
paratively small effect on the standard deviation of range error o^ryg? especially in the 
two-star measurement method. In fact, for this method, could be doubled without 

significantly changing the value of crryg- 
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Figure 14.- Accuracy characteristics of various range-measurement methods. 1-sigma angular-measurement error is 10 arc-seconds. 
(Methods involving star-to-body measurements also include angular measurement between Earth and Moon.) 


In figure 14 the curve for the angular -diameter measurement method shows this 
method to be accurate only near the Earth. The other three curves in the figure apply 
to the determination of range by either equation (13) or by the perturbation method 
(eq. (14) or (15)). In the perturbation method an additional error, which is predictable 
under certain conditions, is caused by out-of-plane effects and is discussed in the next 
section. The rapid rise indicated near the Earth for these curves is caused by the close 
proximity of the vehicle to the Earth-Moon line as explained in reference 9. The impor- 
tant items to note from figure 14 are (1) that the range -determination error (at first mid- 
course) is somewhat lower for star-to-Moon measurements than for star-to-Earth 
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measurements and (2) that the error is almost twice as large for two-star measurements 
(or stable-platform measurements) as for one-star measurements. 

The star-to-body angle would be most likely determined from measurements such 
as those shown in figure 15. 


Earth 



Figure 15.- Measurements for 
determining O, 

Under the assumption that the Moon is more nearly spherical than the Earth, error due 
to oblateness effects would be much smaller if the Moon were used as the body. 

Out-of-Plane Effects on Range Measurement 

As previously stated, the perturbation method of measuring range has an inherent 
error due to out-of-plane effects. This error is caused by the fact that for any perturbed 
trajectory the vehicle will not lie in the nominal Earth-Moon-vehicle plane at Tpf. For 
the two-star measurement method, as well as for the one-star measurement method in 
which the star is in or very near the nominal plane, this error is essentially predictable 
as shown in figure 16. The plots in figure 16 were obtained from the data presented in 
table VI. (For calculations involving stars much out of the nominal Earth-Moon-vehicle 
plane, an additional factor is required in the last term of eqs. (18) and (19) as shown in 
reference 9. This factor has been included in the calculations for table VI but is impor- 
tant only for the out-of-plane star Alpheratz (6 = 28.391°).) In figure 16 the ordinate is 
the difference between the true value of Ar^g and the value calculated from equa- 
tion (14). Values of Arye - Arye,calc shown in table VI for various stars, most of 
which are fictitious, to show how the error is affected by star location. 

The empirical curves in figure 16 can be used to correct the range measurement. 
The scatter in figure 16(b) is less than that in figure 16(a), the implication being that 
the use of a star in the direction of the Moon-Earth line (or the two-star method) 
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TABLE VI.- ERRORS IN RANGE DETERMINATION BY PERTURBATION METHOD 


No measurement error 


Nv Condition 

^^ve> 

Perturbed-^\^ km 

trajectory ' 

no. 

Acye ~ ^^ve^calcj “ 

Star in nominal 
Earth-Moon- 
vehicle plane; 

e = 90®, 

6 = 0° 

Star in nominal 
Earth-Moon- 
vehicle plane; 
e = 64.242®, 

6 = 0® 

Star in direction of 
Earth-Moon line; 

6 = 20.957®, 

6 = 0® 

(*) 

Star 10 arc-tnin out 
of Earth-Moon- 
vehicle plane; 
e = 64.2419®, 

6 = 0.16670 

Star 20 arc-min out 
of Earth-Moon- 
vehicle plane; 

9 = 64.2422®, 

6 = 0.3333® 

Star 1.0® out of 
Earth-Moon- 
vehicle plane; 
0 = 64.246°, 

6 = 1.0® 

Star in direction of 
major axis of 
error ellipsoid; 

9 = 39,036®, 

6 = 4.698® 

Regulus; 

9 = 67.751®, 
6 = 6.922® 

Alpheratz; 
6 = 67.523®, 
6 = 28.391® 

1 

135.87 

0.43 

0.43 

0.55 

-0.07 

-0.07 

-3.08 

-18.15 

-23.42 

-108.34 

2 

271.85 

2.48 

2.60 

2.74 

1.34 

0.21 

-4,55 

-34.68 

-45.22 

-214.85 

3 

-342.54 

1.50 

1.62 

1.62 

1.50 

1.37 

-0.92 

-9.85 

-2.90 

-18.68 

4 

-686.23 

6.93 

7.07 

6.93 

6.80 

6.55 

5.64 

-15.74 

-1.89 

-32.79 

6 

365.09 

-0.44 

0.06 

1.58 

-3.47 

-7.01 

-24,09 

-135.41 

-143.06 

-653.05 

8 

-840.68 

11.35 

11.49 

11.84 

13.24 

15.01 

22.15 

65.08 

84.08 

344,45 

9 

836.58 

11.51 

13.04 

12.13 

9.99 

8.22 

1.09 

-41.64 

-60.18 

-315.26 

10 

-419.37 

2.72 

3.84 

2.86 

3.73 

4.49 

8.06 

29.42 

38.94 

168,51 

20 

-696.20 

4.60 

4.71 

4.85 

5.99 

7.13 

12.30 

45.01 

56.52 

242.05 

21 

-631.79 

3.50 

3.63 

3.74 

4.76 

5.90 

10.37 

39.95 

49.90 

215.55 

22 

-299.68 

-0.41 

-0,29 

-0.41 

-0.29 

-0.16 

-0.02 

1.69 

2.29 

11.34 

23 

-963.57 

14.71 

14.71 

14.84 

15.72 

16.74 

20.93 

36.31 

57.11 

209.95 

24 

1049.16 

20.41 

20.65 

21.43 

18.13 

15.73 

5.52 

-56.13 

-81.76 

-444.17 

25 

-1055.57 

20.45 

20.70 

21.46 

23.23 

25.63 

35.86 

97.96 

124.06 

495.55 

26 

1184.33 

28.63 

29.00 

30.26 

25.85 

22.81 

10.23 

-66.32 

-97.04 

-542.96 

27 

-1192.06 

28.51 

28,89 

30.28 

31.93 

34.96 

47.56 

124.79 

156.21 

614.50 

28 

-166.44 

0.64 

0.75 

0.90 

-0.38 

-1.51 

-5.67 

-47,93 

-42.43 

-195.31 

29 

-509.71 

5.30 

5.31 

5.42 

4.16 

2.89 

-1.72 

-54.35 

-42.13 

-209.76 

30 

507.40 

5.37 

5.37 

5.49 

6.63 

7.77 

12.39 

66.24 

52.95 

225.01 

31 

-1102.85 

17.06 

17.57 

19.47 

21.24 

24.90 

39.85 

144,37 

168.64 

710.48 

32 

1095.55 

1 17.47 

17.85 

19.74 

14.18 

10.64 

-4.29 

-108.10 

-132.15 

-664.77 

33 

-1286.84 

20.81 

21.82 

25.86 

27.26 

32.69 

54.61 

216.93 

244.33 

879.80 

34 

1276.88 

20.53 

21.54 

25.71 

16.10 

10.79 

-11.11 

-171.83 

-200.17 

-986,10 

36 

-592.59 0.09 

0.48 

1.74 

3.51 

6.41 

18.71 

117.43 

123.99 

567.18 

37 

591.47 0.29 

0.67 

1.94 

-2.36 

-5.40 

-17.57 

-115.35 

-123,02 

-562.80 

38 

-1263.70 25.64 

25.76 

26.15 

27.54 

29.31 

36.83 

73.27 

100.97 

372,36 

39 

-1609.62 41.47 

41.60 

41.98 

43.37 

45.01 

52.08 

78.13 

112.69 

371.14 

40 

1591.37 41.07 

41.19 

41.45 

39.55 

37.91 

30.87 

6.32 

-39.31 

-268.89 


These results also apply to two-star method (stars not necessarily in the Earth-Moon-vehicle plane). 



would lead to the most accurate range determination. In the most accurate determina- 
tion the maximum deviation from the empirical curve is about 5 kilometers, which adds 
little (in a statistical sense) to errors with magnitudes of those shown in figure 14. It 
can be seen from table VI, however, that any star within 10 to 20 arc-minutes of the 
nominal Earth-Moon-vehicle plane would suffice for range measurement. For the star 
10 arc-minutes out of this nominal plane, the maximum spread about the empirical curve 
would be about 7 kilometers, except for trajectory 34 which would be 10 kilometers. For 
the star 20 arc-minutes out of the plane, the maximum spread would be about 10 kilo- 
meters, except for trajectory 34 which would be 15 kilometers. The errors become 
much more pronounced for stars more than 20 arc-minutes out of the plane and cannot 
be accounted for empirically. 

Summary of Range-Determination Accuracy 

The important facts concerning manual range determination at first midcourse are: 

1. If the one-star method is employed, the perturbation procedure must be used 
because only the deviation (from the nominal) in the angle can be measured. 

2. Even though it is inherently more inaccurate than the one-star method, the two- 
star method will probably be required for most lunar missions since preflight analysis 
would show no navigation stars close enough to the nominal plane to facilitate use of the 
one-star method. 

3. In the two- star method the total angle can be measured, and thus range can be 
calculated directly from equation (13). 

SUMMARY OF CALCULATIONS REQUIRED IN GUIDANCE PROCEDURE 

The preflight and onboard calculations required in the present guidance procedure, 
along with the pertinent equations, are briefly summarized in the following outline. No 
attempt is made herein to systematize the procedure for the easiest, fastest, and most 
accurate calculation method. Furthermore, the outline does not include the equations or 
steps required for alining the thrust vector for the guidance maneuver; they are simple 
to derive (for example, see ref. 12). 

The prefUght calculations are made for the nominal trajectory and apply to the 
predetermined nominal time of the measurements Tpf. For the preflight calculations 
it is assumed that the nominal values of the position and velocity coordinates with respect 
to both the Earth and the Moon, as well as the Earth-Moon distance and the direction 
cosines of the required stars, are known and that the required transition matrices have 
been calculated. 
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Required preflight calculations: 

1. Angle A cos A = 


Tye + ^vm ^em 




2. Angle B (or C) sin B = 


^vm 

i^eni 


sin A 


Note: The nominal values of A and of B (or C) are required only if the 
perturbation method of measuring range is used. 

3. Star-to-body angles { 6 i; 62 ; O 3 ) cos 9 p- 

4. Derivatives of A, 9\, 62 , and 03 

with respect to time for updating 

the measurements Equations (Bl) and (B3) in reference 9 


5. Slopes from plots similar to figure 7 SE^ETpfJ ^^/^Tpf 


6 . Guidance matrices 


[e]: [g]; [h]) 


Equations (7) and ( 8 ) 


Equation (21) in reference 9 


Required onboard calculations: 

1. Convert all measurements to 

a common time Tpf 

2. Range rye Equation (13), (14), (15), (18), or (19) 

3. Values AA, AB, AC, or A0, if the perturbation 

method to calculate range rye is used For example, AA = A^ - An 

4. Quantities D, E, and F for Tpf Equation (1) 

5. Multiply quantities D, E, and F by Kj, K 2 , and K 3 , respectively 

6 . Midcourse-velocity corrections AV and AVf 

(AVf can be calculated at any time during the 

mission prior to aim-point time^ Equations (9) and (10) 


GUIDANCE ACCURACY CHARACTERISTICS 


In this section the more important errors associated with the guidance procedure 
are defined and analyzed. It is assumed that there are no errors due to human limita- 
tions in calculating numbers with many significant figures. Although the guidance pro- 
cedure described herein has been termed "Manual," all necessary onboard calculations 
do not have to be done by longhand or slide rule. Longhand calculation would require too 
much time, and there is no conceivable guidance procedvme in which sufficient accuracy 
could be produced by calculations made entirely by slide rule. Slide-rule accuracy is 
sufficient for range (or position) determination; but, for the matrix multiplications 
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required in the guidance computations, at least six-place accuracy is necessary (for a 
guidance procedure with measurement errors of magnitudes shown in fig. 14). Hence, it 
is assumed that some means of performing simple desk-type calculations would be avail- 
able to the navigator. 


Test Conditions 

In order to determine the general accuracy of the guidance method, specific condi- 
tions were selected to be representative of most lunar missions. 

First mid course maneuver .- The time of the first midcourse maneuver Tj was 
selected to be 10 hours. This time is close to that for optimum fuel requirements, 
although this proximity is not necessary for making the subsequent error analysis. The 
time of the position fix Tpf was selected to be 9.5 hours, and thus one-half hour was 
allowed for performing the guidance calculations (eq. (9)) and for orienting the thrust 
vector. For practical purposes this time interval seems reasonable; however, the 
ensuir^ results could also apply to a longer (or shorter) time interval. A trial calcula- 
tion of a first-midcourse -maneuver velocity correction was performed in approximately 
15 minutes; no checking procedures were included. For this calculation it was assumed 
that the various measurements required for the position fix had been made shortly before 
Tpf and had all been corrected and updated to Tpf. (A simple method for updating 
these measurements is given in ref. 9.) The calculation for the second-midcourse - 
maneuver velocity correction (eq. (10)) could be made any time during flight prior to the 
final (aim-point) time Tf. 

Aim points .- In the error analysis two different aim points are investigated: one 
at perilune (Tf = 70 hours, 37 minutes) and one at a point near the lunar sphere of influ- 
ence (Tf = 56 hours). The variation of the errors between these two points is also 
included. For the aim point at perilune, the first midcourse maneuver is used to correct 
the actual (perturbed) trajectory to the nominal perilune (position only). The second 
midcourse maneuver (made at perilune) is required to correct the inherent velocity error 
due to the first midcourse maneuver. For the aim point near the lunar sphere of influ- 
ence, the actual trajectory is corrected to the point on the nominal trajectory 56 hours 
from nominal injection time, where the second midcourse maneuver is applied. 

The elements of the transition matrices depart from linearity as the central body 
(Moon) is approached. (See refs. 6 and 13, for example.) Hence, for comparison pur- 
poses, the two aim points were investigated, with the one at 56 hours chosen far enough 
away from the Moon to be relatively unaffected by it. There are other reasons for use of 
different aim points. For instance, for a particular type of mission, it may be simpler 
and/or more accurate to apply a terminal guidance procedure (to a trajectory that is 
supposedly on the nominal course) after having made the second-midcourse -maneuver 
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correction prior to reaching the Moon (or planet), whereas, for another type of mission, 
the second-midcourse-maneuver correction might be easily included in the terminal 
guidance maneuver or in the deboost maneuver. 

Other conditions for which to aim, such as a given lunar approach asymptote, could 
be incorporated in the present method. This procedure would involve a set of terminal 
coordinates which would linearize the transition matrices for a wide range of 
perturbations. 


Transition Matrices 

Examples .- The transition matrices required for the guidance calculations are 
shown in table VII. These matrices apply to the nominal trajectory (see table n) and 
were generated by the program in reference 7. The 3x3 submatrices in the first, sec- 
ond, and last matrices ((a), (b), and (e) parts of table vn) are used in the guidance equa- 
tions. The fourth matrix ((d) part of table VTI) is the transformation matrix required for 
mapping the covariance matrix of injection errors (T = 0) to Tpf = 9.5 hours (required 
in error-ellipsoid calculations previously discussed). The third matrix ((c) part of 
table VII) was not used in the present analysis, but is included for possible use if an 
error analysis between the two times shown is desired. 

Linear approximation effects .- In order to obtain a general idea of the error caused 
by linearization of the transition matrices, "perfect" guidance (i.e., guidance with no 
measurement or execution errors) was applied to the 40 perturbed trajectories in 
table ni. (The fuel requirements for these trajectories are summarized in a subsequent 
figure.) The guidance was performed for the two different aim points and the results are 
given in the last four columns of the table. These results apply to the nominal trajectory 
of table II which comes to within 1665 kilometers of the Moon's surface 
(rp = 3403 kilometers). As expected, the discrepancies at the aim points are larger for 
the trajectories that are farther from the nominal trajectory. The curves in figure 17 
illustrate the progressive increase in the aim-point errors as the trajectory is perturbed 
farther from the nominal. These curves are shown for the position and velocity com- 
ponents for which perturbations have the largest effect on the trajectory (for example, 
compare the first 12 trajectories in table HI). 

As noted from table HI, the errors at the perilune aim point are larger than those 
at the 56-hour aim point, with the larger percentage increase being in the velocity error. 
The position error increases up to three times whereas the velocity error increases as 
much as 100 times. The great increase in the errors at perilune is caused largely by 
the aforementioned effect of the Moon on the assumption of linearity. All position and 
velocity errors shown in table HI, however, are generally small in comparison with those 
caused by errors in the measurements required for the guidance calculations. (For 
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TABLE Vn.- STATE-TRANSITION MATRICES 


jjlxponent given at end of each numberj 


(a) Matrix for transitioning from 10 to 56 hr 


0.18397068 00 

0.16367418 00 

0.90780576 -01 

0.14164047 +06 

0.25834212 

+04 

0,13015172 

+04 

0.20618665 00 

0.24776638 +01 

0,12994105 +01 

0.27734888 +04 

0.20458872 

+06 

0.36354531 

+05 

0.11394173 00 

0.12931567 +01 

0.87041364 00 

0.14036259 +04 

0.36264995 

+05 

0.15882821 

+06 

0.60015734 -05 

0.19635148 -05 

0.10805260 -05 

0.75236952 00 

0.86392488 

-01 

0.42457492 

-01 

0.17664494 -05 

0.14366116 -04 

0.11632002 -04 

0.68260383 -01 

0.15076761 

+01 

0.47507985 

00 

0.95564674 -06 

0.11469357 -04 

-0.13010716 -06 

0.32313430 -01 

0,47209395 

00 

0.89892318 

00 


(b) Matrix for transitioning from 

10 hr to perilune (70 

hr, 37 min) 




-0.16804177 00 

0,11607129 +01 

0.63233023 00 

0.37760592 +06 

0,53896308 

+05 

0.37585866 

+05 

0.26125418 00 

0.22270348 +01 

0.13125725 +01 

-0.48636477 +05 

0.19298204 

+06 

0.56686290 

+05 

0.12585353 00 

0.12528484 +01 

0,58871565 00 

-0.85384851 +04 

0.57612417 

+05 

0.10891148 

+06 

-0.76817188 -04 

-0.47622401 -03 

-0.27299891 -03 

0.35488797 +02 

-0.40694930 

+02 

-0.16977675 

+02 

0.53477561 -04 

-0,66364833 -04 

-0.45173942 -04 

-0.79632870 +02 

-0.35366604 

+01 

0.89227447 

+01 

0.12195098 -04 

-0.13885297 -03 

-0.41910182 -04 

-0.31335628 +02 

0.76927085 

+01 

-0.29492002 

+02 


(c) Matrix for transitioning from 

56 hr to perilune (70 

hr, 37 min) 




0.21204965 +01 

0.16041629 00 

0.13682995 00 

0.10168238 +06 

0.24281123 

+03 

0.31475286 

+04 

-0.28738580 00 

0.65319876 00 

0.80186812 +01 

-0.16681536 +05 

0.37706714 

+05 

0.37513230 

+04 

-0.54127618 -01 

0.97726759 -01 

0.48629285 00 

-0.40828978 +04 

0.44341775 

+04 

0.29211317 

+05 

0.20071121 -03 

-0.13350981 -03 

-0.42088363 -04 

0.10793003 +02 

-0.80042390 

+01 

-0.26208864 

+01 

-0.44804018 -03 

-0.30181057 -04 

0.47614757 -04 

-0.21763512 +02 

0.14839254 

+01 

0.36260822 

+01 

-0.17590183 -03 

0.59466203 -04 

-0.16495501 -03 

-0.84971374 +01 

0.41743676 

+01 

-0.76179589 

+01 



(d) Matrix for transitioning from 0 to 9.5 hr 




0.41948920 +02 

-0.61193071 +02 

-0.30832064 +02 

0.88249186 +05 

0.28259802 

+04 

0.61275881 

+04 

0,16772236 +02 

-0.54822736 +02 

-0.20628127 +02 

0.67458479 +05 

0.20433145 

+05 

0.12345455 

+05 

0.12472478 +02 

-0.24840537 +02 

-0.29174238 +02 

0,41960052 +05 

0.88120858 

+04 

0.11975222 

+05 

0.10350756 -02 

-0.17212393 -02 

-0.88055056 -03 

0.24779977 +01 

0.22689039 

00 

0.25883674 

00 

0.94137598 -03 

-0.24844021 -02 

-0.11336515 -02 

0.31798932 +01 

0.76216386 

00 

0.57231656 

00 

0.59758215 -03 

-0.12763631 -02 

-0.10274083 -02 

0,18953488 +01 

0.41093051 

00 

0.37390006 

00 


(e) Matrix for transitioning from 10 to 9.5 hr ^e., [0]g g 

J H" ■ 

hr 

Ml/ 

) 


0.99962930 00 

0.14864207 -03 

0.83659335 -04 

-0.17997740 +04 

-0.91610628 

-01 

-0.51570603 

-oil 

0.14864130 -03 

0.10004870 +01 

0.49733545 -03 

-0,91610389 -01 

-0.18002968 

+04 

-0.30332312 

00 ' 

0.83658908 -04 

0.49733550 -03 

0.99988375 00 

-0.51570472 -01 

-0.30332313 

00 

-0.17999289 

+04 j 

0.41830011 -06 

-0.16975658 -06 

-0.95562211 -07 

0.99961767 00 

0.15689809 

-03 

0.88340428 

-04 

-0.16975439 -06 

-0.55001961 -06 

-0.56190313 -06 

0.15689728 -03 

0.10005027 

+01 

0.51402240 

-03 

-0.95560998 -07 

-0.56190326 -06 

0.13135906 -06 

0.88339983 -04 

0.51402245 

-03 

0.99987967 

00 
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Off-nominal velocity, u, m/sec 
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(b) Effect of error in injection velocity. 

Figure 17.~ Examples of errors at aim point (Tf = 56 hr) caused by linear 
approximation made in using transition-matrix theory. Results are from 
guidance maneuver at Ti == 10 hours. 



I 


example, see table I.) Hence, for the range of injection-condition perturbations shown in 
table ni (magnitudes considered to be representative of manned injection systems) and 
for the nominal trajectory used herein, errors caused by linearization of the transition 
matrices would have only minor effects on the overall guidance accuracy. 

Fuel Requirements 

As a matter of interest, the fuel requirements (i.e., midcourse velocity corrections 
AV and AVf) for the present method are shown in figure 18. The data shown apply to 
"perfect" guidance maneuvers made for some of the trajectories in table ni and are for 
the aim point at Tf = 56 hours. For the aim point at nominal perilune, the first- 
midcourse velocity corrections are about 2 percent lower than those shown for the aim 
point at Tf = 56 hours and the second-midcourse velocity corrections are from 5 to 
55 percent lower, with the majority being about 40 percent lower. 



Figure 18.- Midcourse-maneuver fuel requirements for the aim point at Tj = 56 hours. 

The data in figure 18 show that for the perturbed trajectories used in the present 
study, the fuel requirements are essentially a direct function of Ar at Tpf = 9.5 hours. 
This relationship occurs because, for these trajectories, Ar ® s where s is the 
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position deviation at Tpf = 9.5 hours. For trajectories with larger injection errors, 

Ar -f s; hence, the direct relationship between fuel requirements and Ar would not hold. 

Effect of Measurement Error 

The results of an analysis to determine the effect of measurement error on guid- 
ance accuracy are shown in figures 19 to 28. Trajectory 10 was close to the nominal 
trajectory and was used to insure that none of the resulting guidance error could be 
attributed to linear approximation effects (transition matrices). 

In the error analysis the effect of scatter in the plots of 6D as a function of D 
(fig. 7) was not included inasmuch as the point for trajectory 10 fell directly on the 
faired line. The point could be off the line for any given perturbed trajectory; however, 
the additional guidance error incurred would be relatively small statistically. 

The three stars simulated for 
the angular measurements in the 
guidance equation (that is, for deter- 
mining D, E, and F in eq. (1)) 
are indicated in figure 19. Other 
combinations of stars could be used 
for these measurements, as long as 
each star produces an adequate plot 
of 6D as a function of D (small 
amount of scatter). If one or two of 
these stars were also to be used for 
determining the range, then each 
must also meet the requirements 
previously discussed under "Range 
Determination. " 

Error at first midcourse. - 
0 10 20 30 ifO 50 Data in figures 20 and 21 depict mea- 

Angle between star and center of Earth SUrement errors at Tnf = 9.5 hOUrS 

at T f = 9.5 hours, 0, deg ^ 

^ and the corresponding errors in the 

Figure 19.- Effective error in range due to measurement error of first-midcourse -maneuver velocity 

10 arc-seconds in star-to-Earth angle at Tp, = 9.5 hours. Tl = 10 hours due to these mea- 

surement errors. 

The measurements required for the guidance procedure, as shown by equation (1), 
are range r^^ and the three star-to-body angular measurements di^H) ^2,a> ^3,a* 

Actually, as previously discussed, the range measurement is composed of a combination 
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of angular measurements, but it will be treated as a different type of measurement for 
error-analysis purposes. The range -measurement error is the predominant error since 
most of the position-determination error is in the direction of range (for example, see 
ref. 9 ); therefore, the guidance accuracies are presented for the range -measurement 
error only, as well as for all the measurement errors. The results pertaining only to 
range -measurement error are shown because they give a good first-order approximation 
of the overall guidance accuracy characteristics. 


In figure 20, the first-midcourse -maneuver velocity errors are represented by 
open symbols for errors in range measurement only (at Tpf = 9.5 hours) and by solid 
symbols for errors in all measurements. The open-symbol data were obtained by 
merely solving the guidance equation (eq. (9)) with the use of discrete values of range 

error. As indicated in the figure, 
the solid-symbol data were obtained 



Figure 20.- First-midcourse-maneuver velocity error at Tj = 10 hours 
due to measurement errors at Tpf = 9.5 hours. Open symbols are 
for range-measurement error only, whereas Monte Carlo results also 


include angular-measurement errors 
10 arc-seconds each. 


Oq , and On of 

03 


from a Monte Carlo procedure. In 
the Monte Carlo operation, 100 
samples of measurement errors 
were used to simulate 100 position- 
determination errors at 
Tpf = 9.5 hours. Each sample con- 
sisted of an error in each of the 
four measurements r, 6 i, 62 , 

and 63 . Each error was assumed 
to be uncorrelated and from a nor- 
mal (random) distribution of errors. 
The covariance matrices obtained 
for the first-midcourse -maneuver 
velocity errors are shown in 
table Vni. Both of the actual 
matrices (see parts (a) and (c) of 
table vm) were determined by per- 
forming the calculation 

< Ay > |ax Ay AzJ 
Az 


for each Monte Carlo sample, summing the respective elements for the 100 samples, and 
then dividing each summation by 100 to get the average value. 
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TABLE Vm.- COVARIANCE MATRICES OF FIRST -MID COURSE -MANEUVER 
VELOaTY ERRORS DUE TO MEASUREMENT ERRORS 


[Exponent given at end of each numbe0 

(a) Actual matrix; Oj. = 10 km, (Jq^ = 10 arc-sec, = 10 arc-sec, = 10 arc-sec 



X 

y 

z 

X 

y 

z 

X 

0 

0 

0 

0 

0 

0 

y 

0 

0 

0 

0 

0 

0 

z 

0 

0 

0 

0 

0 

0 

X 

0 

0 

0 

0.159076 -06 

-0.610913 -07 

-0.121657 -06 

y 

0 

0 

0 

-0.610913 -07 

0.309944 -06 

0.170459 -06 

z 

0 

0 

0 

-0.121657 -06 

0.170459 -06 

0.206078 -06 


(b) Approximated matrix; = 10 km, (Jq^ = 10 arc -sec, a ^2 = arc-sec, 

= 10 arc -sec 



X 

y 

z 

X 

y 

z 

X 

0 

0 

0 

0 

0 

0 

y 

0 

0 

0 

0 

0 

0 

z 

0 

0 

0 

0 

0 

0 

X 

0 

0 

0 

0.22503 -06 

0 

0 

y 

0 

0 

0 

0 

0.22503 -06 

0 

z 

0 

0 

0 

0 

0 

0.22503 -06 


(c) Actual matrix; 

Up = 22 km. 

- 10 arc-sec, 

a ^2 arc-sec. 

= 10 arc-sec 


X 

y 

z 

X 

y 

z 

x 

0 

0 

0 

0 

0 

0 

y 

0 

0 

0 

0 

0 

0 

z 

0 

0 

0 

0 

0 

0 

X 

0 

0 

0 

0.221381 -06 

-0.158528 -06 

-0.206143 -06 

y 

0 

0 

0 

-0.158528 -06 

0.142303 -05 

0.768859 -06 

z 

0 

0 

0 

-0.206143 -06 

0.768859 -06 

0.572972 -06 


(d) Approximated matrix; Op = 22 km, cjq^ - 10 arc-sec, Oq^ = 10 arc-sec, = 10 arc-sec 



X 

y 

z 

X 

y 

z 

X 

0 

0 

0 

0 

0 

0 

y 

0 

0 

0 

0 

0 

0 

z 

0 

0 

0 

0 

0 

0 

X 

0 

0 

0 

0.739128 -06 

0 

0 

y 

0 

0 

0 

0 

0.739128 -06 

0 

z 

0 

0 

0 

0 

0 

0.739128 -06 
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The relative effect of errors in the range measurement and in the three star-to- 
body angular measurements is shown in figure 21. The circular symbols were deter- 
mined from figure 20. For example, 
for CTj. = 10 kilometers, the total 
error in determining the off-nominal 
position components D, E, and F 

(eq. (1)) is (lo2 + 10^ + lO^)^/^ kilo- 
meters. The square symbol in fig- 
ure 21 was determined from fig- 
ure 19; that is, the position- 
determination error was obtained by 
taking the square root of the sum of 
the squares of the three values shown 
in figure 19. The values shown in 
figure 19 are actually errors in D, 

E, and F and were determined from 
equation (1). It should be noted that 

P 20 4o 6o 8o error can be reduced by using 

_ stars for which Q is close to 0*^ 

|_'^D J ^ (or 180^^). Also, the increment in 

midcourse -maneuver velocity error 

Figure 21.- First-midcourse-maneuver velocity error at Tj = 10 hours due to error in the three angular 

due to total error in determining off-nominal position components D, measurements (at cr = 10 kilometers 

E, and F at T,,. = 9.5 hours. ^ 

and Or = 22 kilometers) in figure 20 

is about identical to the value shown 

in figure 21 for angular-measurement 

errors only. 

The important result shown in figures 20 and 21 is the fact that rar^e- measurement 
error is the predominant error. Hence, every practical consideration should be used to 
reduce this error. Methods to accomplish this reduction are discussed subsequently. 

Er ror at aim point .- The trajectory position and velocity aim-point errors due to 
measurement error are shown in figures 22 to 26. Some of these statistical results are 
shown on normal-probability plots in figures 27 to 28. The position errors (figs. 22 
and 24) and perilune-distance errors (fig. 28) are shown for both correlated and uncor- 
related errors. 
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As in figure 20, the results in figures 22 to 25 are shown for range -measurement 
error only as well as for errors in all measurements. As previously noted, error in the 
range measurement is the predominant error and analysis of this error, separately, pro- 
vides a first-order approximation of the overall accuracy at the aim point. The results 
for this condition were easily determined by incorporating the corresponding first- 
midcourse -maneuver velocity errors (discrete) of figure 20 into the perturbed trajectory 
and propagating it to the aim point. The relative effects of position error and velocity 
error at Tf = 56 hours (caused by measurement error) on the position and velocity 
error at the Moon are shown in table IX. The position error at Tf = 56 hours accounts 
for 2/3 of the error at the Moon, probably because of the change in the flight-path angle 
(with respect to the Moon) at the lunar sphere of influence produced by this position 
change. The curves in figures 22 to 25 are faired to zero, the error due to linearization 
being disregarded. As shown by table HI, this error would be relatively small and would 
depend on the particular perturbed trajectory used in the error analysis. 


TABLE IX.- COMPARISON OF THE EFFECTS OF POSITION AND VELOCITY ERRORS 
AT Tf = 56 HOURS ON ACCURACY AT PERILUNE 


Results shown for perturbed trajectory with error of -10 km in range measurement at 




At nominal perilune time (Tf = 70 hr, 37 min) — 

Condition at Tf = 56 hr 

Range from 

Position deviation 

Selenocentric 

Velocity deviation 


Moon center, 

from nominal, s, 

velocity, 

from nominal, u, 


km 

km 

m/sec 

m/sec 

Position and velocity errors* 

3534.6011 

142.785 

1953.8167 

27.2331 

Velocity error only (position error 
assumed to be zero) 

3442.8394 

41.469 

1972.6556 

8.4478 

Position error only (velocity error 
assumed to be zero) 

3494.9175 

101.511 

1961.4204 

19.2872 

Nominal (no errors) 

3403.6367 


1980.5083 



*s = 140.375 km; u = 1.0617 m/sec. 
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0 10 20 30 ho 50 


One -sigma error in range measurement, km 

Figure 22.- Effect of range- measurement error on position accuracy at 
Tf = 56 hours. Monte Carlo results also include angular-measurement 
errors Oo , Oa , and of 10 arc- seconds each. 

®2 


The accuracy characteristics are shown in figures 22 to 25 for the two aim points 
previously discussed and in figure 26 for any point that might be selected within the lunar 
sphere of influence. The Monte Carlo results in figures 22 to 25, as well as the results 
in figure 26, were obtained by mapping the covariance matrices of first-midcourse- 
maneuver velocity errors (parts (a) and (c) of table VTII) to the .aim points. The mapping 
was accomplished with the use of the program in reference 7. In several instances, the 
100 trajectories resulting from the Monte Carlo samples were individually propagated to 
the aim point and the resulting 1 -sigma error calculated. The results obtained in this 
manner matched those obtained from mapping the covariance matrices, the indication 
being that the sample size of 100 Monte Carlo selections was large enough for a true 
statistical representation of the errors. The results were also obtained by mapping 
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spherical covariance matrices of first-midcourse -maneuver velocity errors (parts (b) 
and (d) of table VTH). These matrices were arbitrarily obtained by taking each diagonal 
element equal to 1/3 the total Oy2 (from fig. 20 or from the sum of the diagonal 
elements of the actual covariance matrices in table Vni) and setting the off-diagonal 
elements equal to zero. Figures 22 to 25 show that use of this simple procedure may be 
adequate when range -measurement error is as small as about 10 kilometers but not 
when it is larger. If the spherical covariance matrix could be used, the 100 Monte Carlo 
selections would not be needed for determining the covariance matrix; the data of fig- 
ure 21 would be sufficient. 

A Monte Carlo error analysis was also made for the method which employs sepa- 
rate measurements at Tpf and Tj. This method differs from the present one in that 
the knowledge of position deviation at Tj is not based on the measurement of the posi- 
tion deviation at Tpf; hence, the position errors at these two times are not correlated. 
The results of the analysis are shown in table I along with some results obtained with 


O 

0) 



o 



Figure 23.- Effect of range- measurement error on velocity accuracy at 
Tf = 56 hours. Monte Carlo results also include angular-measurement 
errors , Op , and of 10 arc-seconds each. 

^2 ^3 
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the present method. The data shown correspond to 1-sigma measurement errors of 
10 kilometers in the range and 10 arc-seconds in the star-to-body angles. Based on the 
comparison in table I, the present method is about 15 times as accurate as the two- 
star measurement method. 

It may be noted from figure 24 that the results at perilune are the same whether 
the aim point is at 56 hours or at the nominal perilune time (70 hours, 37 minutes). 

This fact is true for both the Monte Carlo results and the results of using errors only in 
the range measurement in figure 24, for all the velocity results in figure 25, and for the 
perilune-distance-error results in figure 28. Hence, the selection of the aim point has 
no effect on the final accuracy at the Moon. Also shown in figures 22 and 24 for the 
position error is the fact that the effect of error in the three star-to-body measurements 
is reduced as Or increases. 



One -sigma error in range measurement, cr^, km 


Figure 24- Effect of range- measurement error on position accuracy at 
perilune. Monte Carlo results also Include angular-measurement 
errors and Oq^ of 10 arc-seconds each. Results apply 

to Tf = 56 hours and = nominal perilune time. 
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Figure 25.- Effect of range- measurement error on velocity accuracy at 
perilune. Monte Carlo results also include angular-measurement 
errors Oq , On , and of 10 arc-seconds each. Results apply 

C/2 

to Tf = 56 hours and = nominal perilune time. 


In comparing the accuracies at the two aim points it is seen that there is little 
difference in the position values, whereas the velocity error at perilune is much larger 
than that at Tf = 56 hours. The large errors at perilune, however, are not too critical 
since they could be easily corrected for in a terminal-guidance procedure. In figure 26, 
the bends in the curves for times near perilune are caused by the focusing effect of the 
Moon. 

The position-error results from the Monte Carlo analysis for Tf = 56 hours are 
presented on a normal-probability plot in figure 27. Similarly, perilune-distance-error 
results are plotted in figure 28; these results apply to either of the aim points. Normal- 
probability graph paper can be used as a simple qualitative test for a normal distribu- 
tion; the scale is adjusted so as to produce a straight-line probability plot for normally 
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Figure 26.- Aim- point errors showing focusing effect of 
Moon on trajectory error resulting from midcourse- 
guidance uncorrelated measurement error. Data also 
Include angular-measurement errors and o^ 

of 10 arc-seconds each. 
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Figure 27.- Probability distribution of the position devia- 
tion at Tj = 56 hours for 100 Monte Carlo samples of 
correlatea position-measurement errors at first mid- 
course. Or = 10 km; Og^ = 10 arc-sec; 

Og^ = 10 arc-sec; Og^ = 10 arc-sec. 


position errors at the Moon are not 


distributed data. Examination of the plots 
in figures 27 and 28 shows that the aim-point 
errors can be considered normally distrib- 
uted; therefore, they can be analyzed by 
standard statistical procedures. For 
example, approximately 68 percent of the 
errors will be less than la and the maximum 
error expected would be about 3a. 

The probability values in figure 27 
were obtained by assuming an equal distri- 
bution of "negative" position deviations at 
the aim point. The position-deviation value 
indicated for la in this figure compares with 
that shown for the corresponding case in 
figure 22. It is of interest to note from 
figure 28 that the 1-sigma values of the 
perilime -distance error are considerably 
different from the values of position devia- 
tion shown for corresponding cases in fig- 
ure 24. This difference indicates that the 
rily in a radial (from the Moon) direction. 


Effect of correlated erro rs.- As seen by equation (1), if one of the stars used for 
^1> ^ 2 ) ^3 used for determining the value of ra, then there is 100 percent 

correlation between the errors in ra and 0 ^ g^, &2 a’ ^3 a> ^ different stars are 

used, there is no correlation. Depending upon where the star is located, the correla- 
tion can have either a helpful or detrimental effect on the accuracy of D, E, or F. 

The effect will be detrimental if the location of the star is such that the two factors in 
ra cos 0a both increase (or decrease) for a given error in the measurement of 0a- 
For example, if 0j a error by +10 arc-seconds, then cos 0j a decreased; and 

if, for example, 0j a ^® used in equation (19) to determine range, the value for ra 
will also be decreased for certain locations of the star. 


In the error analysis (see figs. 22 to 28), the results are shown for uncorrelated 
and correlated errors. For the correlated errors, the star was located so that the 
effect of the correlation was detrimental; that is, the position deviation at the aim points 
increased by 18 percent (figs. 22 and 24) and the per ilune -distance error, by about 
15 percent (fig. 28). Effects of correlated errors that are equal and opposite to those 
shown in these figures could be obtained for other stars. 
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Figure 28,- Probability distribution of perilune-distance error for 100 Monte 
Carlo samples of position-measurement errors at first midcourse. Results 
apply to Tf = 56 hours and Tj = nominal perilune time, o^ = 10 km; 

O 0 J = 10 arc-sec; = 10 arc-sec; og^ = 10 arc-sec. 

Effect of range -measurement error .- The range-measurement error that can be 
expected for 1-sigma errors of 10 arc-seconds in the angular measurements is about 
40 kilometers and 80 kilometers for methods using one and two star measurements, 
respectively. (See fig. 14.) As shown in figures 22 to 25, range -measurement error is 
the predominant factor in determining the guidance accuracy and errors of this order 
would produce excessive errors at the aim point. Any means to reduce the range- 
measurement error would be of great benefit. One possibility would be to Incorporate 
into the measurement schedule a series of range measurements. More specifically, 
each angular measurement required for range determination would be measured a num-, 
ber of times (updated to a common time) and averaged so that the resulting accuracy in 
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range determination would be increased by the statistical law 

_ ®^unave raged 
'^averaged ^ 

where N is the number of measurements made for each angle. 

There is a limit as to how far the reduction in range-measurement error should 
be taken. This limit is set by the inherent error due to scatter in the plots of 6 q as a 
function of D (fig. 7), which is equivalent to a measurement error. From the equation 

D = Ta cos 0a ” ^n ®n = ^a <^0® ®a ~ constant 

it is seen that 

dD = cos 0a dr - ra sin 0a d0 

For 0a near 0° (or 180°) and for representative errors in r and 0 of 10 kilometers 
and 10 arc-seconds, respectively, it can be shown that the error in ra is approximately 
equal to the error in D. It is seen in figure 7 that plots with minimum scatter have a 
scatter of at least ±1 kilometer in 6D and that this scatter represents an error of about 
±10 kilometers in D or in ra. It is not possible to reduce the minimum scatter shown 
in figure 7 by use of other stars; therefore, any attempt to reduce the error in ra much 
below 10 kilometers by adding more measurements would be useless. 

Effect of Nominal Perilune Altitude 

The results in this report pertain to trajectories with a perilune altitude of 
1665 kilometers — that is, trajectories designed to miss the surface of the Moon by a 
relatively large distance of 1665 kilometers. There may be some question as to what 
extent the guidance accuracy would be affected for trajectories aimed closer to the Moon. 
Any changes in the accuracy would be caused by (1) effect of linear approximation in 
transition-matrix theory and (2) Moon’s focusing effect on error propagation. In order 
to resolve this question, the accuracy characteristics were investigated for a nominal 
trajectory similar to the one used throughout this report, but with a perilune altitude of 
155 kilometers. 

The effect of linear approximation on aim-point error was examined for trajecto- 
ries with perturbations corresponding to those for the seven trajectories from table HI 
listed in table X. The results shown in table X for the two nominal perilune altitudes 
indicate that this effect does not increase appreciably for trajectories closer to the Moon. 
The effect on the velocity error, however, is seen to increase in most of the trajectories 
by a factor of 2 or 3 . 
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TABLE X.- EFFECT ON ACCURACY OF LINEAR APPROXIMATION USED IN 
TRANSITION-MATRIX THEORY FOR TWO NOMINAL 
PERILUNE ALTITUDES 



Off-nominal conditions at nominal perilune time for perfect 


measurements and perfect guidance maneuver at first 

Trajectory with same 


midcourse when — 


perturbation as 





trajectory no. 

Nominal perilune altitude = 1665 km 

Nominal perilune altitude = 155 km 


s, km 

u, m/sec 

s, km 

u, m/sec 

4 

45.028 

10.676 

36.211 

18.818 

8 

40.700 

9.589 

57.819 

31.126 

10 

13.804 

3.005 

41.109 

22.279 

26 

29.669 

7.922 

27.782 

18.684 

29 

39.704 

9.070 

44.648 

23.414 

33 

82.143 

20.390 

105.453 

60.190 

36 

36.697 

9.009 

37.779 

21.314 


The Moon's focus factor (see appen- 
dix C) is shown in figure 29 as a function 
of nominal perilune altitude. This factor 
causes the so-called Moon's focusing 
effect on the position and velocity errors 
very near the Moon, as illustrated in fig- 
ure 26. The focusing effect on the errors 
for the near-miss nominal trajectory was 
investigated by mapping the actual covari- 
ance matrices of first-midcourse - 
maneuver velocity error (parts (a) and (c) 
of table Vni) to the nominal perilune time. 
The aim-point errors for this trajectory 
are shown in table XI along with those for 
the nominal trajectory having the larger 
perilune altitude. The difference in the errors for. the two altitudes is due, for the most 
part, to the focusing effect. As shown, this effect on the position errors is about the 
same for both altitudes; on the other hand, the effect on the velocity errors is much 
greater for the lower perilune altitude. 



Figure 29.- Effect of perilune altitude on focus factor of Moon. 
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TABLE XI.- PERILUNE AIM-POINT ERRORS FOR TWO NOMINAL TRAJECTORIES 



Maneuvering Errors 

The effect of guidance -maneuvering errors on the accuracy at the aim point was 

briefly examined. (See figs. 30 to 32.) 
No statistical analysis was performed; 
rather, the individual effects of the 



0 12 3^5 

Timing error, minutes late 


Figure 30.- Effect of first-midcourse-maneuver timing errors on 
position and veiocity accuracy at Tf = 56 hours. 


three most important errors associ- 
ated with the first midcourse maneuver 
(other than alining the reference 
system) were determined. 

Timing errors .- The effects of 
first-midcourse -maneuver timing 
errors on the position and velocity 
errors at the aim point (Tf = 56 hours) 
are shown in figure 30 for two per- 
turbed trajectories. As shown in 
table ni, trajectory 24 is perturbed 
to a much higher degree than trajec- 
tory 22. The position and velocity 
deviations in figure 30, as well as in 
figures 31 and 32, are given in terms 
of the perfect guidance aim-point con- 
ditions of the actual (perturbed) trajec- 
tory rather than the aim-point condi- 
tions for the nominal trajectory in 
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Deviation in velocity u' at T_ = 5^ hours 





order to eliminate the linear approximation effects. It is apparent by comparing the data 
of figure 30 with those of figures 22 and 23 that, regardless of the amount the actual tra- 
jectory is off the nominal trajectory, a timing error, even on the order of several min- 
utes, has a negligible effect in the overall guidance accuracy. 

Pointing errors .- The effect of pointing-direction errors is shown in figure 31. 

The data were calculated for only one perturbed trajectory (trajectory 24 in table HI) but 
are representative of trajectories which require a vector AV in the general direction 
of the trajectory. For perturbed trajectories requiring a vector AV more nearly per- 
pendicular to the trajectory, the effects of pointing errors would be smaller than those 
shown. The results are given for pointing errors that are either in-plane or out-of-plane 
and indicate that pointing errors up to about 0.5° could be tolerated — that is, would have 
no large effect on the overall aim-point accuracy (based on the errors shown in figs. 22 
and 23 for Oj. = 10 kilometers). 
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Figure 32.- Effect of errors in magnitude of first- midcourse- maneuver velocity 
correction on position and velocity accuracy at Tf = 56 hours for trajectory 24. 
Required first- midcourse- maneuver velocity correction AV is 66.63652 m/sec. 
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Magnitude errors .- The effect of first-midcourse-maneuver velocity-magnitude 
(thrust-cutoff) errors is shown in figure 32. The results are for trajectory 24 of table in 
which is perturbed relatively far from the nominal trajectory. Naturally, for trajectories 
perturbed to a lesser degree, the effect of the magnitude error would be less. Also, the 
results shown are representative of perturbed trajectories which require a velocity cor- 
rection in the general direction of the trajectory. Calculations pertaining to the various 
perturbed trajectories of table ni showed that the direction of the required velocity cor- 
rection was almost always in the direction of the trajectory. For the trajectory repre- 
sented in figure 32 and for a comparison similar to that made in figure 31, it is evident 
that errors in the first-midcourse -maneuver velocity magnitude up to about 1 percent 
could be tolerated in the present manual guidance procedure. 

CONCLUSIONS 

A method for determination of midcourse -velocity corrections for a space vehicle 
has been presented. The method has been applied to a lunar trajectory; however, it 
would also apply to an interplanetary trajectory. The method makes use of knowledge 
(preflight calculations) of a nominal or reference trajectory and a linear perturbation 
procedure (transition matrices) to determine the velocity corrections. The linearity 
characteristics of the method are such that sufficient accuracy can be obtained even 
though large differences may exist between the actual and nominal trajectories. In addi- 
tion to an instrument to measure angles, the only items required by the navigator are a 
small amount of precalculated information on the nominal trajectory, several simple 
equations which can be solved within a time of 15 minutes by employing a desk-type cal- 
culator, and some manual procedure to aline the thrust vector. 

It has been shown that the present method essentially involves a two-step calcula- 
tion procedure: (1) calculation of the range to the Earth by use of one or two star-to- 
body measurements along with the angle measurement between the Earth and Moon, and 
(2) calculation of the midcourse-velocity vector by use of the range determination and 
three star-to-Earth measurements. From the error analyses made on range determina- 
tion and on the aim-point conditions, the important results concerning manual guidance 
are; 

1. Error in the range measurement is the predominant factor in the overall guidance 
accuracy. The aim-point accuracy is essentially proportional to this error, regardless 
of any other measurement error or maneuver error. 

2. For a 1-sigma error of 10 arc-seconds in angular measurements, the error in 
determining range for one set of required measurements is about 40 and 80 kilometers, 
respectively, for the one-star and two-star methods. Errors of this magnitude lead to 
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relatively large aim-point errors; however, the accuracy can be increased greatly through 
use of several range measurements. 

3. Because of the inherently large range -measurement error and its effect on the 
aim-point accuracy, more work on improving range -measurement accuracy is needed. 

4. The aim-point position error does not change greatly as the Moon is approached, 
whereas the aim-point velocity error increases rapidly near the Moon. This increased 
error, however, can be easily corrected for in a terminal guidance procedure. 

5. The final accuracy at the Moon (nominal perilune time) is independent of the 
aim point selected in the guidance procedure. 

6. The effect of correlated errors caused by making measurements to certain stars 
can be advantageously used to increase the guidance accuracy. 

7. The aim-point errors have been shown to be normally distributed; hence, stan- 
dard statistical procedures can be applied in analysis of these errors. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., June 6, 1967, 

125-17-05-09-23. 
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APPENDIX A 


A METHOD OF SELECTING INJECTION ERRORS 

The errors in both the position vector and the velocity vector at injection are 
basically of two types: (1) direction errors caused by misorientation of the vector and 
(2) magnitude errors. For example, in considering the velocity errors (see fig. 33), the 
components of the injection velocity are 

X = V cos S' cos (j)' 

y = V cos 6' sin <p' i (Al) 

z = V sin 6' 

J 


Z 



Figure 33.- Orientation of injection-velocity vector. 


After taking the variation in equations (Al), the velocity error is 

6x = ^ 6V - z cos 69' - y 6(p' 

6y = ^ 6V - z sin 0' 60’ + x 60’ > 

6z = — 6V + V cos 0’ 60’ 


(A2) 
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APPENDIX A 


Similar equations would be obtained for error in the injection-position vector. In solving 
equations (A2), if 6V, 69', and 6(j)' are assvuned to be independent random errors, the 
variances of the errors in the velocity components are 





J 


The values used for Oy, and in the present analysis were obtained as follows. 

A value of 10 m/sec was arbitrarily chosen for Oy. Then, in order to determine values 
for (Tg, and this value of 10 m/sec was applied separately to each of the three 

nominal velocity components and the largest difference in each angle (0* and (^') from 
the nominal value was selected (see table XU). 


TABLE Xn.- DATA USED IN DETERMINING Og, AND FOR EQUATIONS (A3) 


Condition 

■ 

X, 

km/sec 

y, 

km/sec 

z, 

km/sec 

V, 

(x2 ^ y2 ^ 

km/sec 

sin 6', 
[VV| 

deg 

sin <p\ 
|y/V cos O'j 

0’, 

deg 

Nominal 


-10.519851 

-2.6309587 

-2.0283715 

11.031930 

0.1838637 

10.594890 

0.2426221 

14.041349 

Error applied to 

X 

-10.509851 

-2.6309587 

-2.0283715 

11.022395 

.1840228 

10.604164 

.2428393 

14.054178 

Error applied to 

y 

-10.519851 

-2.6209587 

-2.0283715 

11.029550 

.1839034 

10.597199 

.2417531 

13.990029 

Error applied to 

z 

-10.519851 

-2.6309587 

-2.0183715 

11.030096 

.1829877 

10.543831 

.2426221 

14.041349 


As shown in table xn, the maximum differences in 9' and <p' from their nominal 
values are, in both instances, approximately 0.051°. By using this value (in radians) for 
Oq, and in equations (A3), the values, 

0jj2 = 0.00009948 (km/sec)^ 

= 0.00009355 (km/sec)^ 

a^2 _ 0.00009653 (km/sec) ^ 

as listed in table IV, are obtained. 
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APPENDIX B 


TWO-STAR METHOD OF RANGE DETERMINATION 

As described in the section "Range Determination," a measurement is required of 
the angle B or the angle C (depending on the method used) or of the deviation of this 
angle from the nominal value. Once this angle is obtained, procedures for determining 
range are relatively simple. The angle B (or the angle C) can be obtained from 
onboard measurements of two star-to-body angles. 


To star 1 



Figure 34.- cones of position formed by star-to-body measurements. 
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Figure 34 represents the two cones of position formed by the two star-to-body 
measurements. The vehicle will lie on one of the two lines of position (intersection). 

The geometry shown in figure 34 (that is, stars in a direction from the vehicle opposite 
to that of the Earth (or Moon)) represents star-to-body angles between 90° and 180°. For 
angles less than 90° the stars would be in the same general direction from the vehicle as 
the Earth (or Moon). The angle B (or C), as shown in figure 12, can be determined 
from the calculation of the direction of the pertinent line of position, along with informa- 
tion on the direction of the line between the Earth and the Moon. Although the calcula- 
tions are lengthy, they can be performed before the mission and presented in graphical 
form. 

The derivation of the direction of the line of position follows and is given in two 
parts. First, the coordinates of the line are determined for a right-hand rectangular 
coordinate system in which the Z"-axis is along the line of sight to one of the stars. 

Next, these coordinates are transformed to the conventional inertial coordinate system 
used throughout the paper. 

In figure 34 the right-hand rectangular coordinate system is centered in a plane 
perpendicular to the axis of cone 1 and at a unit distance from the body center; the 
Z"-axis is along the axis of this cone and the Y"-axis intersects the axis of cone 2. The 
intersection of the X'Tf "-plane with the two cones forms an ellipse and a circle, as shown 
in figure 35. 



Figure 35.- Ellipse and circle formed by intersection of two cones and X"Y"-p!ane. 
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A limitation to this approach is the requirement that 

(/3 + r) < 90^ 

(See fig. 36.) Thus, in selecting stars to be used in this method of range determination, 
the two axes of the cones should be reasonably close together; that is, jS < (90^ - y). 
This condition, in general, is easily met, as the same stars are ordinarily used for the 
position-component measurements in the guidance equations; as shown by table V, the 
best stars for these measurements provide lines of sight which are all fairly close 
together. 

The equation for the circle is 

xft2 _j_ ytt2 _ tan^Q? (Bl) 

and the equation for the ellipse is 

^ + (B2) 

a2 b2 


Solving these equations simultaneously yields the intersection of the two curves, given by 

\l/2 


„ _ -a^c ±(a^b2 - a2b2ta n2o! + c2a2b2 - a2b^ + b^tan^o;) 


2 1,2 
a - b 


(B3) 


X' 




tan'^of - 



(B4) 


As noted, the equations yield four possible combinations of values for x" and y". The 
correct sign of y" is easily determined, however, by testing equation (B4); the wrong 
sign gives no solution. The correct sign of x" must be obtained from predetermined 
knowledge of the direction of the line of position for the nominal trajectory. The line of 
position for the actual trajectory will lie close to this nominal line. As previously 
stipulated, 

z" = 1 (B5) 

The terms a2, b2, and c can be solved by the procedure explained in the discussion 
that follows. 

First, the coordinates in the equation for cone 2 (see fig. 36) are given in terms of 
the X", Y", and Z" system by 
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or 


X"’ 

Vi 

Z»it 

s- J 


1 

0 

0 


0 

cos 0 
-sin /3 


0 

sin ^ 
cos ^ 




yiti _ yit QQQ ^ ^ 2 " sin ^ 
z'" = -y" sin ^ + z** cos ^ 


r-Cone 2: x*"^ + y*"2 = z’"^tan2y 



Figure 36.- Two cones in two rectangular right-hand coordinate systems. 


In the X", Y", and Z" system, then, the cone equations are given as 
Cone 1: x”2 + y"2 - z"2tan2o! - 

Cone 2: x"2 + (y" cos ^ + z" sin 0)^ = (-y" sin 0 + z" cos ^)2tan2y 
After expanding the equation for cone 2 and collecting terms, 
jj)t2 ^ y”2(cos2^ - sin2^ tan2y) + 2y”z” sin 0 cos ^(1 + tan2y) = z"2(cos2/3 tan2y - sm^0) 
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A plane perpendicular to the Z"-axis at a unit distance from the origin is given by the 
equation 

z" = 1 

The intersection of this plane with cone 1 is given by 

x't2 + yit2 _ tan2o! 


and the intersection with cone 2 is given by 

x"2 + y"2(cos2/3 - sin2/3 tan2y) + 2y" sin /3 cos ^(i + tan2y) = cos2/3 tan2y - sin2^3 

Adding the term sin /3 cos ^(1 + t an y) both sides of the equation for the intersec- 
cos2^ - sin2/3 tan^y 

tion of the plane with cone 2 gives 


x”^ + (cos 2|3 - sin^/3 tan2y) 


y„2 ^ 2y"sin /3 cos g(l + tan2y) 
cos 2|3 - sin2/3 tan2y 


sin2)3 cos2g(l + tan2y)^ 


(c6s^/ 3 - sin^/3 tan^y) J 


= cos 


2^ tan2y - sin2^ + sin2^ cos2^(l + tan2y)^ 
cos 2/3 - sin2/3 tan2y 


or, after some algebraic manipulation, 


x”2 + (cos 2/3 - sin2/3 tan2y) 


,,, sin g cos g(l + tan2y) 


-,2 


cos^/3 - sin^jS tan^y 


tan2 


y 


cos^;8 - sin^/3 tan^y 


or 


cos^/3 - sin^^ tan2y ^,,2 , (cos2^ - sin2/3 tan2y)^ 
r -X + 

tan^y 


tan^y 


sin j3 cos )3(1 + tan2y) 


cos^/3 - sin^/3 tan^y 


= 1 


which is in the form of equation (B2). Therefore, 


q 2 = tan2y 

cos 2/3 - sin2/3 tan2y 
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b2= if 

cos^P - sin2jS tan2y 
^ _ sin j3 cos /3(1 + tan2y) 

cos 2(3 - sin2/3 tan^y 

Next, the coordinates given by equations (B3), (B4), and (B5) must be transformed 
to the X, Y, Z axis system. If i, j, and k are unit vectors along the X-, Y-, and 
Z-axes, respectively, then 

z" = Iji + mjj + njk 

where Ij, mj, and nj are known direction cosines of star 1 with respect to these 
axes. 

Since by definition (fig. 37) X" is perpendicular to the plane passing through the 
lines of sight to the two stars, x" is given by the cross product of the position vectors 
of the two stars as 


X" = (min2 - nim2)i + (nil2 - lin2)j + (lini2 - mil2^k 

where I 2 , m 2 , and n 2 are known direction cosines of star 2 with respect to the X-, 
Y-, and Z-axes. 



Figure 37.- Illustration of two rectangular right-hand coordinate systems. 
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With the vectors and defined, the vector is then given by 

yT» sz X 


or 


y" = - l2nii2 _ l2ni^ + linin2^i 

+ ^nin2mi - ni2ni2 - li^m2 + lil2mi^f 
+ - ^ 2 ^ 1 ^ + nimini2^k 


The direction cosines of the vectors x", y", and z" form the transformation 
matrix 0 such that 


/ N 

X 




>=CpJN 

y"> 

bJ 


z'' 
L J 


where the elements of [P] are defined by 


Pll 


min2 - njm2 


j(min2 - nim2)^ + (nil2 - lin2)^ + (lim2 - 


1/2 


“1^2 - 


X' 


Pl 2 


’^1^2 " ^1^2 


X' 


P13 


ll“2 - “1I2 


X' 


P21 


(ll™l“2 - l2“l^ - hH^ + h"in2) 


[7 2 2 

- l2mi2 - l2ni2 + linin2) + (nin2mi - m2m2 - li^ni2 + lll2mi) + (I1I2IU - -"2“!^ + nii"l'"2)J 


lim^mg - l2mj2 . l^nin2 
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P22 = 


nin2ini - + l^lgin^ 


P23 = 


lll2^1 - ^ 1^2 - + njinimg 


P31 = h 


P32 = mi 
P33 = ni 

The angle B (or C) , then, is determined by the dot product of the position vec- 
tors of the vehicle and Moon (or Earth) center, as referenced from the Earth (or Moon) 
For example. 


XeyXem + YevS^em + ^igyZem 
cos B 

^ev^em 

where Xem, Yem* ^em known ephemeris values. 
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FOCUS FACTOR OF MOON 
The focusing effect of the Moon is shown in figure 38. 



Figure 38.- Sketch showing focusing effect of the Moon. 

In this figure, a is the semimajor axis and e is the eccentricity of the orbit. The 
distance b<» is the effective aim-point radius — that is, the radius within which a 
desired perilune distance or radius rp can be reached. It is seen that 

boo = i-p^l + 7“^)‘^os(0oo - D = 

Now, as shown by the well-known equation 

r= I 

1 + e cos 6 
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where 1 is the semilatus rectum, r tends to infinity as 1 + e cos 6 approaches 0, or 
as d approaches the value 



Therefore, from equations (Cl) and (C2), 


(C2) 


boo 



^^p = fFp 


where ff is the focus factor. 

The eccentricity can be determined from 

where p is the product of the universal gravitational constant and mass of the Moon. 
The perilune velocity can be determined from 




where Vg is the vehicle velocity (with respect to the Moon) at the sphere of influence 
and rg is the distance or radius of the sphere of influence (approximately 57 500 km). 

In calculating the results shown in figure 29, data given in reference 14 were used 
to determine the value of Vg for the nominal trajectory used in this report. All results 
shown in the figure apply to trajectories having the same energy (that is, Vg is 
constant). 
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